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ABSTRACT 

We have investigated the color-magnitude diagram of tu Centauri and find that the blue main 
sequence (bMS) can be reproduced only by models that have a of helium abundance in the range 
Y = 0.35-0.40. To explain the faint subgiant branch of the reddest stars ("MS-a/RG-a" sequence), 
isochrones for the observed metallicity ([Fe/H] —0.7) appear to require both a high age (~ 13 
Gyr) and enhanced CNO abundances ([CNO/Fe] ~ 0.9). Y « 0.35 must also be assumed in order to 
counteract the effects of high CNO on turnoff colors, and thereby to obtain a good fit to the relatively 
blue turnoff of this stellar population. This suggest a short chemical evolution period of time (< 1 
Gyr) for u> Cen. Our intermediate-mass (super-)AGB models are able to reproduce the high helium 
abundances, along with [N/Fc] ~ 2 and substantial O depletions if uncertainties in the treatment of 
convection are fully taken into account. These abundance features distinguish the bMS stars from the 
dominant [Fe/H] rj —1.7 population. The most massive super-AGB stellar models (Mzams > 6.8 M Q , 
■^He.corc > 1.245 M ) predict too large N-enhancements, which limits their role in contributing to 
the extreme populations. In order to address the observed central concentration of stars with He-rich 
abundance we show here quantitatively that highly He- and N-enriched AGB ejecta have particularly 
efficient cooling properties. Based on these results and on the reconstruction of the orbit of u> Ccn 
with respect to the Milky Way we propose the galactic plane passage gas purging scenario for the 
chemical evolution of this cluster. The bMS population formed shortly after the purging of most of 
the cluster gas as a result of the passage of uj Cen through the Galactic disk (which occurs today 
every ~ 40Myrs for ui Cen) when the initial-mass function of the dominant population had "burned" 
through most of the Type II supernova mass range. AGB stars would eject most of their masses into 
the gas-depleted cluster through low- velocity winds that sink to the cluster core due to their favorable 
cooling properties and form the bMS population. In our discussion we follow our model through 
four passage events, which could explain not only some key properties of the bMS, but also of the 
MS-a/RGB-a and the s-enriched stars. 

Subject headings: globular clusters: individual (u> Cen) — Hertzsprung-Russell diagram — stars: AGB 
and post-AGB — stars: abundances — stars: evolution — stars: interiors 



1. INTRODUCTION 

Omega Centauri provides an especially valuable con- 
straint on our understanding of chemical evolution be- 
cause of its unusual properties and its proximity, which 
enables us to study its component stellar populations 
in great detail over a very extended range in luminos- 
ity. While it has been known for some time that the 
giants in this syste m encompass a range in [Fe/H] from 

2.2 to -0. 5 (e.g. |Brown fc Wallerstein|1993|[Suntzeff 

& Kraft 1996), the extensive spectroscopic surveys car- 
ried out in the past decade, in particular flSmith et al. 
20001 |HlIke"r et al.||2004| |Kayser et al.||2006F van Loon 



et~aTl|20071 1 Johnson et al.||2008| | Johnson fc Pilachowski 
2010| jMarino etHT. 2012 and references therein), have 
established that the metallicity distribution rises sharply 
from [Fe/H] « -2.2 to a strong peak at [Fe/H] w -1.7, 
with a long tail that drops off to higher metal abun- 
dances containing secondary peaks at [Fe/H] w —1-4, 
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— 1.1, and —0.7. Type la supernovae appear to have 
contributed to the chemical makeup of only the most 
metal-rich stars given that their measured [cc/Fe] abun- 
dance ratios are significantly reduced from the constant 
value of ss 0.35 found in stars haying [Fe/H] < —0.8 
( |Pancino et"al]|2002[ [Origlia et al.|)2003[ ). The elevated 
a-element abundances over most of the range in metal- 
licity (also see Kayser et al.) implies that Type II su- 
pernovae were the major producers of these elements. 
Moreover, the rapid rise in the [La/Fe] ratio between the 
most metal-deficient and the [Fe/H] > — 1.5 populations 
(by a factor of 3 according to Johnson & Pilachowski) 
indicates that s-processing, and therefore, intermediate- 
mass asymptotic-giant branch (AGB) stars contributed 
significantly to the chemistry of uj Cen after the domi- 
nant [Fe/H] = —1.7 population had formed. 

Accompanying these spectroscopic advances have been 
equally impressive improvements in the photometric 
data. The detailed, and deep, co lor-mag nitude dia- 



grams (CMDs) derived by|Lee et al 



Wallersteinl (I2000D ; |Rey et aTT fl2004 ) ; Bcd fnet al. | f2004|) ; 
al. ( 



(|1999); [H^ g5eTfc| 



TTTlanova et al. 
lini ct al. 



Sollima et al. (2005a 
Calamida et al ] ( |2009p 

among others, have also established the existence oFsev- 
eral discrete stellar populations in uj Cen. The most baf- 
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fling one of them is a so-called "blue main sequence" 
(bMS), discovered by (|Anderson||1997|), that is clearly 



separated from, and bluer at a given magnitude than, the 
MS associated with the dominant [Fe/H] ss —1.7 popu- 
lation (see the CMDs reported by, e.g., Bedin et al. and 
Villanova et al.). What has made this di s covery so puz- 
zling is the determination by Piotto et al. ( 2005 ) that the 
former is more metal rich than the latter by S 0.3 dex. 
As discuss ed by P i otto e t al., and anticipated by Bedin 
et al. and Norris ( 2004[ ), the most obvious (only?) way 
to reconcile these observations with stellar evolutionary 
theory is to infer that the bMS stars have unusually high 
helium abundances (Y > 0.35). 

Although it has not yet been possible to make a defini- 
tive connection of the bMS thro ugh the subgiant region 
of the CMP to the RGB ( but se elKing et al.|2012| ) , |John 



son & Pilachowski (120101) have found that the giants in 



their sample with extreme" abundances (i.e., those with 
low C and O, together with high N, Na, and Al; see 
their Fig. 23) constitute a similar proportion of all gi- 
ants (« 27%) as the fraction of all MS stars that are lo- 
cated on the bMS. Furthermore, since (for the most part) 
their radial distributions are quite similar, it is tempting 
to conclude that the aforementioned giants are the de- 
scendants of stars that had occupied the bMS earlier in 
their evolutionary history. However, because the O-poor 
stars do not exhibit a correlation with [La/Fe], while the 
O-rich stars do, and because there is no apparent corre- 
lation of the La abundance with radius, in contrast to 
the observed trend for the O-poor stars, Johnson & Pi- 
lachowski suggest that the observed light element abun- 
dances reflect both (i) the retention by u> Cen of the 
ejecta of AGB stars and (ii) in situ mixing on the RGB. 

There is a second feature of the CMD that is poten- 
tially quite challenging to explain, and that is the exten- 
sion of the reddest giant branch (independently discov- 
ered by |Lee et al.| ( |1999[ ) and by jPancino et at] ( |2000[ ), 
who gave it the name "RGB- a" that has si nce been used 
to refer to it), through the subgiant branch dFerraro et al 



2004 



2010 



to very faint magnitudes on the MS (|Bellini et al. 



Because this is is the most metal-rich of the dis- 
crete populations that have been identified, it is presum- 
ably also the youngest one. Consequently, its age pro- 
vides a key constraint on the timescale over which most of 
the chemical evolution took place in w Cen. Because the 
turnoff of this fiducial sequence is so faint, the age of this 
stellar population must be quite old, implying that the 
bulk of the star formation occurred over a rather short 
interval of time. Indeed, the reason why some investi- 



gations (e.g. Hilker et al.||2004 Rey et al 



2004) derived 



an extended star formation history (~ 3-4 Gyr) in their 
analyses is that they had mistakenly adopted a bright 
turnoff for the RGB-a population. However, it remains 
unclear whether the different stellar populatio ns of ui Cen 



have the same age to within a few xl O 8 yr (D'Antona 
et al. 2011a Valcarce fc Ca tclan 20 ill), or t hey span a 



range m age of ^> 2 Gyr ( Sollima et al H2005a|) or as muc h 
as 3-5 Gyr ( |Stanford et ai1|2006| I Villanova et aL]|2007[ ). 

Among stellar evolution models of different mass 
ranges, intermediate-mass and super-AGB stars have 
been proposed as the source of the He-enriched and oth- 
erwise extreme abundance patterns associated with the 
bMS. However, models of these types of stars do not 
always reproduce the required high He, high N, and 



low O abundance associated with the bMS in a quan- 
titative wa y. For example, the intermed iate-mass AGB 
models by Ventura & D'Antona (2009 specifically the 
Z = 6 x 10 -4 case that is applicable to the first, most 
metal-deficient generation in w Cen) predict Y = 0.329 
and [N/Fe]= 1.84 for M ZA ms = 5.0 M and Y = 0.360 
and [N/Fe] = 1.36 for M ZAM s = 6 . M^ . The super-AGB 



models by Ventura & D'Antona (2011 ) produce a maxi- 
mum ejecta He abundance of Y = 0.36, but information 
for N is not available. The envelope He abundance at the 
end of the s econd dredge-up in the super-AGB models by 
|Siess| ( [2007 ) reach 0.37 for non-overshoot ing models, and 



0.38 for models with core overshooting. Siess (2010) do 
not provide a model for Z — 6 x 10~ 4 , but interpolating 
between their super-AGB thermal pulse calculations for 
Z = 0.001 and Z = 10~ 4 yield He mass fractions in the 
ejecta of Y ~ 0.33, while N is enhanced by a factor of 
10 to 100. Moreover, those models with the largest in- 
creases in the nitrogen abundance also have the smallest 
O depletion, and models with O depletions of ~ 1 dex 
have N enhancements of only ~ 1 dex. While existing 
models clearly point qualitatively in the right direction, 
it is not clear to us what it would take to reproduce the 
"extreme" abundance mixtures that are associated with 
the bMS in uj Cen in a quantitative sense. 

Producing He-abundances in intermediate-mass AGB 
ejecta that exceeds the He-abundance of the "extreme" 
abundance mix in any significant way appears to be very 
difficult within model calculations. This does not leave 
a lot of room (if any) for dilution of the AGB ejecta be- 
fore the formation of the blue main-sequence stars. Dilu- 
tion has been suggested, however , in order to reproduc e 
the abundance anti-correlations (D'Ercole et al. 2008). 
We are considering here rather the origin and evolution 
of one sub-population at a time. The presentation of 
Na-0 abundances by sub-popula tions identified with a 
four-criterion cluster analysis by Gratton et al. (2011 
Fig. 6) rather suggests that the anti-correlation is trie re- 
sult of the superposition of discrete Na and O abundance 
in dif feren t sub-populations, and as we will argue later 
(Sect. 5.2) dilution may not be needed. We therefore in- 
vestigate the evolution of first-generation intermediate- 
mass AGB and super-AGB stars with the goal to gen- 
erate the "extreme" abundance mix without taking into 
account any dilution. 

In any case, the AGB ejecta may preferentially con- 
verge in the cluster c enter via AGB cooling flows 
(D'Antona et al. 2011b and refs. there). Such gas may 
have to be isolated, for example through supernova- 
induced clearing of intracluster gas that is ejected from 
stars having initial masses outside the range that en- 
compasses intermediate-m ass AGB and super-AGB stars 



(Conroy & Spergel 
old stars increases t 



|2011 ), while the tidal stripping of 
tie ratio of second- to first-generation 
stars (Bckki 2011). However, complete SN gas purging 
makes it difficult to envisage how a spread in [Fe/H], like 
that observed in to Cen, is produced. 

The present study has been undertaken to address 
some of the issues described above. In §2, newly com- 
puted sets of isochrones and zero-age horizontal branch 
(ZAHB) loci for different values of Y, Z, and heavy- 
element mixtures are compared with the CMD of uj Cen 
in order to illustrate how the interpretation of the ob- 
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Table 1 

Chemical Abundance Ratios 



[m/Fe] 



m 


St cLiidcird 


Extreme 


lu Ccn a 


1^ V ' V_, ZOUO 


C 


0.0 


-0.7 


-0.5 


-0.7 


N 


0.0 


+2.0 


+1.5 


+2.0 





+0.5 


-0.5 


-0.5 




Nc 


+0.3 


0.0 






Ka 


0.0 


+0.8 


+0.6 


+0.8 


Mg 


+0.3 


0.0 




+0.1 


Al 


0.0 


+ 1.1 


+1.1 


+1.1 


Si 


+0.4 


+0.4 


+0.4 




P 


0.0 


0.0 






S 


+0.3 


+0.3 






CI 


0.0 


0.0 






Ar 


+0.3 


+0.3 






K 


0.0 


0.0 






Ca 


+0.3 


+0.3 


+0.3 


+0.3 


Ti 


+0.3 


+0.3 


+0.3 




Cr 


0.0 


0.0 






Mn 


0.0 


0.0 






Ni 


0.0 


0.0 







a For O-poor giants (jJolmsonfc^ilachows kipOlO^ . 
b For a blue MS star I Bragaglia et al.|2010^ ~ 



servations would be affected by the assumed variations 
in the chemical abundances, and to make some assess- 
ment of the age and helium content of the stars that 
belong to the bMS and MS-a components. §3 investi- 
gates under which assumptions the predicted yields from 
> 5M Q AGB stars can be made to agree with the ob- 
served/inferred high Y, high-N, low-0 abundances. Fi- 
nally, a short summary of the main results of this study 
is given in §4, which also discusses in some detail how 
our stellar evolutionary results together with the efficient 
cooling properties of the ejecta from intermediate-mass 
AGB stars may naturally explain the formation of ad- 
ditional populations like the bMS. A key point in our 
proposed scenario is the periodic purging of the gas from 
globular clusters (or dwarf galaxies) that is expected to 
occur throughout their evolutionary histories whenever 
their orbits cause them to pass through the Galactic 
plane. Implications of this scenario are briefly summa- 
rized in §5. 

2. COLOR-MAGNITUDE DIAGRAM CONSIDERATIONS 

Several grids of stellar evolutionary tracks were com- 
puted using a significantly updated version of the Vic- 
toria code (VandenBerg et al. 2012) which now treats 
the diffusion ot helium (but not the metals) as well 
as extra mixing below envelope convection zones (when 
the y exist) using methods very similar to those described 
by |Proffitt fe Michaud| dl991| ). All of the model com- 
putations reported in this study assumed a value of 
Qmlt = 2.007 for the usual mixing-length parameter, 
so as to satisfy the solar constraint w hen the mix of 
heavy elements derived for the Sun by |Asplund et al.| 
(2009) is assumed. The abundances of the individual 
a-elements in this mixture were increased by approxi- 
mately the amounts found in ve ry metal-deficient stars 
according to Cayrel et al.| ( |2004 [), resulting in the "stan- 
dard" metals mix (see the second column of Table[l]) that 
applies to normal Population II stars, including (presum- 
ably) the dominant [Fe/H] sa —1.7 component of uj Cen. 



The third column of Table[T]lists the "extreme" heavy- 
element mixture that has been assumed in some of the 
comparisons of observations with models to be discussed 
shortly. Remarkably, multip le main sequences have also 
been discovered in NGC 2808 |Piotto et al |(j2007|), and the 
chemical abundances derived bylBragagha et al. (20101 



for one of its bluest MS stars appear to be quite sum 
lar to those found in what arc believed to be the red 
giant descendants of bMS stars in uj Cen (compare the 
entries in the 4th and 5th columns). To maximize the 
impact of these abundance anomalies on the predicted 
properties of stellar models, we have chosen to adopt the 
largest of the [m/Fe] values (in an absolute sense) that 
have been determined for w Cen and NGC 2808 in our 
"extreme" heavy-el ement mix (column 3) . However, we 



note that, although Marino et al 
the [O /Fe] value varies from —0.6 



2011 ) have found that 
o +0.9 at metallicities 
appropriate to bMS stars, the [CNO/Fe] r atio spans the 
relati vely small range of « +0.4-0.7 (see Marino et al. 



20121) . This suggests that the CNO abundances listed in 
the third column of Table 1 are somewhat more extreme 
than those found i n u> Cen stars having [Fe/H] « —1.5. 
On the other hand, [Marino et aT ( 2012j see their Fig. 3) 



find that the [CNO/FeJ value ranges between +0.7-0.9 
at [Fe/H] ~ —1.0, and as discussed below, such high val- 
ues appear to be necessary to explain the faint subgiant 
branch of the MS-a/RG-a population. 

Also worth mentioning is the fact that opacity data for 
bo th mixtures were gen erated using the code described 
by Ferguson et al.| ([2005) for the low-temperature regime, 
while complemet ary high-temperature ta bles (similar to 
those reported by Iglesias fe Rogers|1996 ) were obtained 
via the OPAL websrt^] Note, as well, that the inter- 
polation program developed by P. Bergbusch, which i s 
similar to that described by VandenBerg et alT] (|2006[), 
but with th e significant improvements (see | VandenBerg) 



et al. 2012), was used to produce all of the isochrones 
considered in this investigation. 

To begin our analysis, we first compare, in Fig.fTl the 
main-sequence segments of several isochrones with the 



lower MS photometry of u> Cen given by Sollima et al. 
d2007ah, based on VLT FORS1 observations. The lat- 



ter have been converted to th e \(B — j?) , Mal-plan e 



M) 



assuming (m 

E(B - V) = 0.136 ( |Schlegel ef^T 
4.07 E(B 



13.70 ( Bellazzini et all 2004 ) 
"199 8J), and 



V) (McCall 2004). 

F 



V), A R = 2.UE(B 
Bellazzini et al. have argued that the |Lub| (j20t)2| de- 
tcrmination of E{B — V) = 0.11 ± 0.01 mag is the cur- 
rent best estimate of the foreground reddening, and that 
may well be the case. However, the main purpose of 
Fig. 1 is to determine the value of Y that is needed to 
reproduce the location of the bMS stars relative to the 
main sequence of the dominant [Fe/H] ss —1.7 popu- 
lation. One could obviously adopt the lower value of 
E(B — V) and then increase the predicted B — R colors 
by 0.026 mag to obtain an identical fit of the models to 
the observations in a differential senseH In fact, fainter 
than Mb — 6.5, the synthetic colors had to be adjusted 

5 See http://opalopacity.llnl.gov 

6 To transpose the models from the theoretical to the ob- 
served plane, we have used the color-T B ff relations derived by 
L. Casagrande (see | VandenBerg et al.|2010} from synthetic spectra 
based on the latest MAKUS model atmospheres | |Gustafsson et aT| 
|2008"L 
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cj Cen 
E(B-V) = 0.136 
(m-M) = 13.70 




13 Gyr isochrones 
[Fe/H] Y Z-Mix 

1.70 0.25 standard 

1.40 0.35 standard 

1.40 0.40 standard 

1.40 0.35 extreme 



0.8 



1 1.2 

(B-R)n 



1.4 



Figure 1. Comparison of the lower main-sequence extensions of 
isochrones for the indicated age and chemical abundances with 
the photometry of a; Cen reported by Sollima ct al. (2007aj). The 
adopted values of E(B — V) and the true distance modulus are 
noted in the upper right-hand corner. 

by S(B - R) = 0.02(M B - 6.5) in order for the [Fe/H] 
= —1.70, Y = 0.25 isochrone (the solid curve) to accu- 
rately reproduce the observed MS over the entire mag- 
nitude range that has been ploteed. (Such a color off- 
set could arise if, among other possibilities, the MARCS 
model atmospheres predict too much flux at short wave- 
lengths or there are systematic errors in the model T e ff 
scale due to the treatment of convection or the surface 
boundary conditions.) 

When such a color correction is applied (to all of the 
model loci), one finds that the bMS stars in cj Cen are 
bracketed by isochrones for Y — 0.35 and Y — 0.40 
(the long-dashed curves) o n the assumption o f the ob- 



-1.4, 



served [Fe/H] value ( 
"standard" Pop. II metals mixture 



Piotto et al.||2005[ ) and the 
We note that Y = 



0.39±0.02 has been obtained by ( |King et aL]2012| ) from a 
fit of isochrones to their improved CMU.) Indeed, these 
results are fully consistent with the findings of P iotto 



et al. (see their Fig. 7) and Sollima et al. ( 2007a| their 
Fig. 9), who performed similar comparisons using com- 
pletely independent theoretical computations. (While we 
have opted to plot 13 Gyr isochrones, mainly because 
their extensions to brighter magnitudes are shown in a 
subsequent figure, the location of the lower MS is ob- 
viously independent of the assumed age.) Interestingly, 
fainter than Mb ~ 5.7, isochrones for the "standard" 
and the "extreme" heavy-element mixtures are essen- 
tially identical when the same values of Y and [Fe/H] 



1 - 




ZAHB loci 

[Fe/H] = -1.4. -1.7 (Y = 0.25) 

[Fe/H] = -1.4 (Y = 0.30, 0.35) 

[Fe/H] = -1.4 (Y = 0.35, extreme mix) 

V) = 0.136, (m-M) = 13.70 

J I I I I I I I I I L_ 



-0.3 0.2 

( m F435W _ m F 



0.4 



0.6 



Figure 2. Comparison of several ZAHB loci for the indicated 
propertie s with the hor izontal-branch stars in uj Cen having Mg < 
3.2 (from [Bellini e t al. 2009J on the assumption of the reddening 
and distance modulus that were assumed in the previous figure, as 
noted. As mentioned in the text, the predicted colors had to be 
increased by +0.04 mag in order for the models to reproduce the 
HB population in the blue tail satisfactorily. 

are assumed (compare the short-dashed curve with the 
long-dashed curve for the same helium abundance). (Of 
course, differences would be evident if photometric filters 
are used that are sensitive to the ab undances of CNO 

Bellini et al.|2010 



and/or other heavy elements; see, e 



Sbordone et al.||2011| |Milone et al.| 2012p As others, in - 



cluding those mentioned above and e.g., Norris| ( |2004[ ), 
have concluded, it does not appear to be possible to 
explain the location of the bMS stars in to Cen (and 
NGC 2808) without invoking high Y — provided that, 
indeed, chemical abundance differences arc primarily re- 
sponsible for the CMD anomalies. Indeed, we explored 
the impact of varying the ab undances of O, Mr, and S i 
using isochrones presented by VandenBerg et al. (2012), 
which allow for variations in the abundances of several 
metals in turn, but we were unable to find a satisfac- 
tory alternative explanation (i.e., other than high Y) for 
the location of the bMS relative to the main-sequence 
fiducial of the dominant [Fe/H] ~ — 1.7 population. 

There is one aspect of Fig. [T] that deserves further 
comment. The bMS stars are clearly separated from 
the dominant red MS at 6.7 < Mb < 7.7, but not at 
lower luminosities, where the two populations appear to 
merge. By contrast, the isochrones remain parallel to 
one another over the entire range in Mb that is con- 
sidered. One can speculate that this may be a color- 
T e ff relations effect; i.e., that at sufficiently cool temper- 
atures, stars having the "extreme" metals mixture are 
redder at a fixed T c g and gravity than stars having nor- 
mal Pop. II abundances. To investigate this possibility, 
it would be necessary to compute proper model atmo- 
spheres and synthetic spectra for the two heavy-element 
mixtures and then compute the fluxes in the various fil- 
ter passbands. Additional work along these lines would 
certainly be worthwhile. 

Before turning to an examination of the upper-MS to 
lower-RGB stars in ui Cen, it is instructive to examine 
its HB population. Fig. [2] compares the HB stars that 
have Mp435\v < 3.2 from the published |Bellini et al.| 
(120091) HST data set with ZAHB loci that have been 



calculated for indicated values of [F e/H] and Y, using 
the numerical methods described by VandcnBerg ct al 
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(2000). As indicated, the same reddening and distance 
modulus that were assumed in the previous figure have 
been adopted here, along with Ap43 5w = 4.081 E(B — V ) 



and j4f625W = 2.637 E(B—V) from |Sirianni et al 
However, besides the reddening adjustment 



(2005) 



;he pre- 
dicted colors had to be corrected by +0.04 mag in or- 
der to achieve a satisfactory fit of the ZAHB loci to 
the steeply sloped blue stars. It is not clear why an 
additional zero-point correction is needed, as the same 
ZAHBs appear to provide an excellent match to the 
F606W, F8UW observations of Anderson et all 12008) 



without having to apply any ad hoc color shiftPJ Be 
that as it may, this concern is not of particular impor 
tance for the present analysis, given that we are primarily 
interested in the comparison of predicted and observed 
CMDs in a differential sense. Note that, to accomplish 
the transformation of the models to the observed plane, 
we have used the color-T e g relation s applicable t o AC S 
photometry that were presented by Bedin et al. ( 2005 ) , 
and kindly provided to us by S. Cassisi (2011, priv. com- 
munication) 



As already shown by Sollima et al. (2005b), most of 
the HB stars in u) Cen having (mF435w — "JF625W > 0.0 
appear to be matched quite well by ZAHB loci for [Fc/H] 
values from ~ —1.7 to ~ —1.4 and Y rj 0.25 (the short- 
dashed curves). Although some of the observed stars lie 
between the two solid curves, which represent ZAHBs 
for [Fe/H] = -1.4 and Y = 0.30,0.35, it is not obvi- 
ous whether they truly have higher helium abundances 
or they have simply evolved to their current CMD loca- 
tions from initial structures on ZAHBs for lower values 
of Y. The most noteworthy result of Fig. [2] is the lack of 
any stars brighter or bluer than the ZAHBs for Y = 0.35 
and [Fc/H] = —1.4 — assuming either the "standard" 
or the "extreme" metals mixtures (i.e., the brightest of 
the solid curves and the long-dashed curve, respectively). 
ZAHB loci for lower [Fe/H] values and/or higher helium 
abundances would be even brighter at a given color than 
the latter. Thus, if u Cen does contain stars having 
Y > 0.35, they must evolve to ZAHB locations well down 
the blue tail. This would perhaps not be too surprising 
given that, at the same age and [Fe/H] value, the turnoff 
mass is predicted to be ss 0.12 solar masses less for stars 
having Y — 0.35 than those for Y — 0.25, nearly in- 
dependently of the metals mixture that is assumed. If 
similar, and significant, amounts of mass loss occur on 
the RGB, then stars having high Y should evolve into 
ZAHB structures that are much hotter than those for 
normal helium abundances. This has already been ap- 
preciated by a number of others (e. g. |Busso et al. | |2007| 
D'Antona fc Caloi|2 008 ) , and in fact, |Cassisi et al | (|2009) 
have suggested that the densely populated clump of HB 
stars at niF435w ~ 18.6 (Mp435w ~ 4.4) in the u> Cen 
CMD may be He-rich stars. 

Fig. [3] makes a number of i nteresting compariso ns of 
theoretical isochrones with the Bellini et al. ( 2009 ) ACS 
data for upper-MS, subgiant, and lower- KGB stars of u> 
Cen. (Only a fraction of the total number of observed 
stars have been plotted, for the sake of improved clarity.) 

7 According to A. Sollima (2010, private communication), the 
HB models that were compared with an eariler reduction of the 
same observations in the Sollima et al. (2005) study also required 
some adjustment of the predicted colors. 
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Figure 3. Comparison of 13 Gyr isochrones with the upper main- 
sequence, subgiant, and lower red-giant-branch populations of oj 
Cen. The short-dashed curves assume Y = 0.25, the "standard" 
mix of heavy elements (see Table [TJ|, and [Fe/H] = —1.7, —1.4, 
—1.0, and —0.7, in the direction from left to right. The solid curves 
were computed for [Fe/H] = —1.4 and the same metals mix, but for 

Y = 0.35 and 0.40: the isochrone for the higher helium abundance 
has a bluer turnoff. The long-dashed loci represent isochrones for 

Y = 0.35, the "extreme" mix, and [Fe/H] = —1.4 and —0.7, in the 
direction from left to right. The dot-long-dashed isochrone assumes 
the "extreme" mixture, as well, but a normal helium abundance 
(Y = 0.250). 



Curiously, it was not necessary to apply any correction to 
the isochrone colors, as derived from the transformations 



by Bedin et al. (20051, in order to obtain quite a good 
match to the observed MS stars, which leads us to wonder 
whether the difficulty noted above concerning the colors 
of ZAHB models should be attributed to a small prob- 
lem with these color-T e g relations (for just warm stars) 
or whether our interpretation of the observed HB is in- 
correct. The former explanation would seem to be the 
most likely one given that, if the predicted colors are not 
adjusted to the red, the core He-burning stars of ui Cen 
would lie below, or redward, of all of the ZAHB loci (as- 
suming that the adopted distance modulus is accurate), 
which would be very difficult to explain. 

The short-dashed curves represent 13 Gyr isochrones 
for Y = 0.25, and [Fe/H] = -1.7, -1.4, -1.0, and -0.7 
(in the direction from left to right), on the assumption of 
the "standard" metals mix. (An isochrone for the lowest 
metallicity and an age of 12.5 Gyr would actually pro- 
vide the closest match to the densest distribution of stars 
which define the brightest subgiant branch (SGB), and is 
therefore our best estimate of the age of the most metal- 
poor population of u Cen.) The RGB segments of these 
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isochrones appear to be too red, by ~ 0.05 mag, which 
suggests that either the model T e g scale is somewhat too 
cool and/or the adopted color transformations for low- 
gravity stars are slightly too red. However, in a relative 
sense, the difference in the colors of the predicted RGBs 
for [Fe/H] = —1.7 and —1.0, at a fixed color, is compa- 
rable to the observed width of the main distribution of 
giant stars. 

Isochrones represented by the solid curves and the 
bluest of the long-dashed loci extend to higher lumi- 
nosities the same i sochrones that were p lotted in Fig. [T] 
As emphasized by Soll ima et al.| ( |2005b[ ), among others, 
variations in the helium content have only minor effects 
on the luminosity and morphology of the SGB. This is 
demonstrated by the fact that the two solid isochrones 
(for Y = 0.35 and 0.40) and the short-dashed isochrone 
(for Y = 0.25) — all of which were computed for [Fe/H] 
= —1.4 on the assumption of the "standard" mix of 
heavy elements — are nearly coincident along the sub- 
giant branch. On the other hand, the isochrones for 
Y = 0.35 and [Fe/H] = —1.4 that has been computed 
for the "extreme" metals mix (the aforementioned long- 
dashed isochrone) has a significantly fainter SGB at a 
same age (due to these models having a much higher 
[CNO/Fe] a bundance ratio). 

Although jKing et al.| ( |2012[ ) have found that the SGB 
tentatively identified as the extension of the of the bMS 
can be matched by isochrones for high Y and a normal 
a-enhanced mixture, the sensitivity of subgiant lumi- 
nosities, at a fixed age , to the total abun dance of the 
CNO elements (e.g., see Cassisi et al.|2008 l, clearly com- 
plicates the determination of the ages of ui Cen SGB 
stars if derived from their locations relative to thcorcti- 



cal isochrones for different metallicities (e.g., Hilker et al. 
2004 Stanford et al.|2007[ ) . In order to have robust ages, 
even m a differential sense, it is of critical importance to 
take into account any star-to-star variation s that exist in 
both [ Fe/H] and [CNO/Fe]. According to |Marino et al. 
(2012), the mean [CNO/Fe] ratio varies from « 0.25 at 
[Fe7H] = -1.9 to w 0.75 at [Fe/H] = -1.0 (which is 
very difficult to understand; see their discussion of this 
point), with a dispersion of about 0.2 to 0.4 dex at a 
given metallicity. Such variations of [CNO/Fe] at a fixed 
iron abundanc e will affect predicted age s at the level of 
- 0.6-1.3 Gyr flVandenBerg et al.|[2012] ). 

Finally, Fig. [3] shows that the most metal-rich of the 
discrete stellar populations in u) Cen has a fainter SGB, 
and an appreciably bluer turnoff, than those of a 13 
Gyr isochrone for Y = 0.25 and [Fe/H] = -0.7 ("stan- 
dard" mix) . Surprisingly, the reddest of the long-dashed 
isochrones, for Y = 0.35 and [Fe/H] = —0.7 ("extreme" 
mixture), does reproduce these observations quite well 
— though the inferred value of Y cannot be determined 
very precisely because the upper MS and turnoff of the 
MS-a population is not very well defined in the photom- 
etry that we have usedj^] On the other hand, isochrones 
for the same metallicity and a normal helium content 

8 The fiducial sequence for these stars is clearly defined in the 
CMD reported by Bellini et al. (2010), but those data are presented 
for filter bandpasses that we are unable to model at the present 
time due to the lack of suitable color transformations. Impor- 
tantly, Bellini et al. find a single, well-defined turnoff for the MS-a 
sequence, which justifies our decision to fit isochrones to roughly 
the middle of this distribution of stars. 



(Y = 0.25) are clearly too red and they fail to reproduce 
the slope of the subgiant branch (note the location of the 
dot-long-dashed curve). Isochrones for high Y do not suf- 
fer from such difficulties. Thus, it appears to be a fairly 
robust stellar evolution prediction that the MS-a/RG-a 
stars have high CNO and He abundances. 
T hese results support previous suggestions. In particu- 



lar, Marino et ah] ( 2011 ) noted that "as with bMS stars, 
the colors and the metallicity of the MS- a are consis 



tent with bei ng pop ulated by He-rich stars ( Norris[2004 
Bellini et al. 20101. Moreover, based on observations 
taken m many different ACS bandpasses, the aforemen- 
tioned Bellini et al. paper suggested that this population 
has "peculiar" CNO abundances. Indeed, the predicted 
age of the MS-a/RG-a stars will be less than the age 
of the universe only if they have high CNO abundances 
(assuming that our estimate of the to Cen distance is ac- 
curate). Furthermore, the slope of the subgiant branch, 
which is well defined in Fig. 3, can be reproduced by the 
models only if high Y is also assumed. As far as our 
contention that a high helium abundances is needed to 
match the turnoff color is concerned, we note that (i) the 
T e flf scale of our models agrees well with that derived for 
solar neighborhood stars (having Hipparcos parallaxes) 
over the entire range in [ Fe/H] encompassed by them (see 
VandenBerg et al.||2010[ ), and (ii) isochrones for close to 
the primordial abundance of helium and standard Pop. II 
metal abundances provide a good fit to the [Fe/H] = — 1.7 
population. Stellar models should be quite reliable in a 
differential sense. 

Our CMD analysis suggests that the different stel- 
lar populations in u) Cen have close to the same age, 
and hence that the chemical evolution in this system 
must have occurred over a relatively short period of time 
1 Gyr). (For the MS-a/RG-a stars to be apprecia- 
bly younger than the universe, they they would neeed 
to have a higher [Fe/H] value and/or higher CNO or 
helium abundances than those assumed in the best-fit 
isochrone.) Especially intriguing is the very real possibil- 
ity that the most metal-rich component has a high helium 
abundance. Is there, then, a connection between these 
stars and the bMS? Before offering some suggestions on 
how the chemical evolution in ui Cen may have proceeded 
(in §4) , we will first examine whether the chemical yields 
from AGB models that form out of gas having the "stan- 
dard" mix of heavy elements (column 2 of Table [I]) will 
be similar to the "extreme" mixture (col. 3). 

3. AGB AND SUPER-AGB STELLAR EVOLUTION WITH 
UNCERTAIN CONVECTIVE MIXING PHYSICS 

Our analysis of the CMD of ui Cen suggests that some- 
how out of the first, and dominant, [Fe/H] ss — 1.7 pop- 
ulation, in which the metals have the "standard" a- 
enriched distribution, at least two distinguishable new 
stellar populations have been generated. These are the 
[Fc/H] —1.4 bMS stars and their counterparts on the 
RGB, which appear to have helium abundances between 
Y = 0.35 and 0.40, and the so-called "MS-a/RG-a" pop- 
ulation at an even higher Fe abundance (« —0.7) that 
seems to have a similarly enhanced He abundance as well 
as enhanced C+N+O, judging from the comparison of 
isochrones with the observed photometry. 

In this section we wish to address the question of 
whether current AGB and/or super- AGB models are able 
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Table 2 

Initial abundance distribution 
for the AGB stellar evolution 
calculations. 
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a Corresponds to "standard" in 
Tablc[7J 

b Corresponds to [Fe/H] = -1.7 
for the a-enhanced mixture con- 
sidered here. 

to produce ejecta out of which the bMS population, and 
possibly the MS-a population, can form. 

3.1. Model assumptions 

We have calculated several intermediate-mass AGB 
and super-AGB stellar evolutionary sequences with an 
initial abundance distribution (see Tables]) that closely 
resembles the abundances found in the (dominant gen- 
eration of stars at [Fe/H] w —1.7. The tabulated mass- 
fraction abundances for the metals are, in fact, equivalent 
to the abundances that have been labelled as the "stan- 
dard " mix in Tableffl and t he adopted isotopic ratios are 
from | Asplund et alj ( [2009] ) . 
We use the MESA stellar evolution code (rev. 2941, 



Paxton et al.||2011| ). The MESA paper by Paxton 
et al. already contains verification cases for a wide range 
of stellar evolution cases, including thermal pulses and 
dredge-up in AGB stars. In addition, we have now also 
run massive AGB models with the same initial abun- 
dances and similar enough physics assumptions as those 
chosen for the grid of models including massive AGB 



stars by Herwig ( 2004b ) obtained with the EVOL stellar 
evolution code, and we find the MESA results again to 
be in reasonable agreement. 

We adopt a customized nuclear network (identified as 
sagb_NeNa.net), which is based on the MESA network 
agb.net. To be more specific, we consider 23 chemical 
species (including the major isotopes) for the elements 
from H to Na, ending with 24 Mg, as well as the pp- 
chain, the CNO cycles, the NeNa cycle, He-burning, the 
13 C(a,n) 16 reaction, and the 12 C + 12 C -^ 20 Ne + a re- 
action. The nuclear reacti on rates have been taken from 
the NACRE compilation ( |Angulo et al|l999[ ). 

For the mixing-length parameter we adopt «mlt = 
1.73, which is obtained from a standard solar model. 



MESA is executed in hydrodynamic mode and sufficient 
artificial viscosity is assumed in order to damp out the 
large velocities that are otherwise predicted to occur near 
the surface during the advanced stages of AGB thermal- 
pulse (TP) evolution (especially during the dredge-up 
phase). In addition to the default mesh-point controls, 
we refine the mesh using the abundances of p, 4 He, 13 C 
and 14 N so that the chemical abundance profiles are al- 
ways well resolved. Additional criteria have been im- 
plemented in order to improve the resolution of He-shell 
flashes and the advance of the thin H-burning shell dur- 
ing the interpulse phase as a function of time. We use 
the atmosphere option ' simple_photosphere ' , and we 
assume a mass-loss rate of logM (in M^/yr) —4.8 
to —4.3, which is obtained using the iBlockerl (119951) 
mass-loss rate that has been incorporated in MESA with 
?j = 5x 10~ 4 . (Although some sequences were followed 
for up to 100 TPs, none of the model stars had lost signif- 
icant amounts of mass by the time the calculations were 
terminated.) We also adopt the OPAL opacities that in - 
clude C- and O-enhanced tables ( Iglesias fc Rogers|1996[ ). 
In the pre- AGB evolution we consider convective bound- 
ary mi xing (CBM) according to the exponential CBM 



mode l (Freytag et al. 1996 Herwig et al. 1997 Herwig 
20001, with. / = U. 014 at all convective boundaries 

We have chosen initial masses of 5.0 M and 7.0 M 
to represent, in turn, an intermediate-mass AGB star 
having a CO core and a hydrogen-free core mass of 
M H = 0.991 M at the first TP, and a super-AGB star 
with an ONe core and a hydrogen-free core mass of 
M H = 1.2677 M Q at the first TP. We present below the 
results that have been obtained for these two cases when 
both standard and modified assumptions are made about 
stellar interior mixing processes. 

3.2. Intermediate-mass AGB star model 
3.2.1. Standard mixing assumption 

We will first describe the surface abundance evolution 
of the standard 5 M Q calculation and how it relates to 
the evolution of convection zones in some detail (Fig.Q. 
For the standard case, we assume convective boundary 
mixing with / co = 0.004 at the bottom of the convec- 
tive envelope and / p d cz = 0.002 at the bottom of the 
pulse-driven convection zone. Both of these values are 
much smaller than those normally assumed in compu- 
tations of low-mass AGB models. As shown by |Her-| 
|wig| ([2005) a CBM parameter as large as indicated by 
s-process-constraints in low-mass AGB stars (/ ce ~ 0.13) 
would lead to very vigorous hot dredge-up which even- 
tually evolves into a corrosive flame penetrating deeper 
and deeper into the core. The model evolution is soon 
terminated in this case due to high mass loss induced by 
extreme luminosities. Such a scenario seems incompati- 
ble with the requirements for the ejecta of intermediate- 
mass and super AGB stars in ui Cen. This issue does, 
however, highlight shortcomings regarding our present 
models of CBM in the deep stellar interior. 

At the beginning of the evolution all chemical species 
have surface abundances that are unmodified from those 
listed in Table[2j and they do not change until the end 
of H-core burning. Then, at log(-fi na i — t) = 7, the first 
dredge-up leads to a small decrease of the surface 12 C 
abundance and a corresponding increase of the surface 
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Figure 4. Surface abundance evolution of dominant CNO iso- 
topes and 4 He of 5Mq stellar evolution tracks with standard mix- 
ing assumptions (top panel) and Kippenhahn diagram (bottom 
panel). The time axis is reverse- logarthmic to accomodate the 
accelerating evolution, ffl na i = 9.077424 X 10 7 yr. 



14 N abundance. No changes are predicted during He-core 
burning, after which major changes, again initially for 
12 C and 14 N but then for 4 He. They are the result of the 
second dredge-up, which starts at log(ifi na i — t) = 5.7 and 
ends at 5.25. The second dredge-up, which takes place 
even before the first thermal pulse on the AGB, mixes 
into the envelope a large fraction of the He that was 
produced by H-shell burning during the He-core burn- 
ing phase. This mixing event is responsible for a ma- 
jor increase in the surface He abundance, from initially 
X( 4 He) = 0.250 to 0.353 for this case. By the end of 
the second dredge-up the 14 N abundance has increased 
by +0.6dex, 12 C has decreased by 0.3dex and 16 has 
decreased by 0.08dex. 

These modifications of the surface, and therefore en- 
velope composition, all point in the right direction, but 
quantitatively they are far from matching the "extreme" 
abundances listed in TablcJT] However, further process- 
ing of the envelope takes place during the thermal-pulse 
AGB phase that starts shortly after the end of the second 
dredge-up. Individual TP-related mixing events, like the 
third dredge-up and the He-shell flash convection zones 
are fully resolved in the stellar evolution calculation, but 



not in the Kippenhahn plot shown in Fig.[4j However, for 
the 5 M case, the repeated action of the third dredge- 
up after each thermal pulse is evident from the periodic 
step-like increase of the 12 C and 16 surface abundance 
evolution, which is also shown in Fig.|4j The thermal- 
pulse AGB evolution starts at log(t final — t) = 5.2, and 
though little or no third dredge-up occurs during the 
five He-shell flashes, there is a steep initial increase of 
the 14 N envelope abundance up to a mass fraction of 
log(A) = —3.0 at log(tfi na i — t) ~ 5.1. This is due, at 
first, to the destruction of 16 O when the third dredge- 
up is weak, and then to just the effects of deep third 
dredge-up (typically Amd up = 1.7 x 10~ 3 M per pulse) 
after each of the next six TPs, These events mix into 
the envelope primary 12 C and 16 from He-shell burn- 
ing. Note that , due to the hot dredge-up in our models 
I Herwig 2004a), the immediate CN-cycling conversion of 
C to N during the dredge-up phase is responsible for 
approximately 20% (a fraction that decreases to below 
10% in subsequent thermal pulses) of the 14 N production 
per pulse cycle (also see Fig. [5] and the discussion below). 
Hot-bottom burning during the interpulse phase is re- 
sponsible for the remainder of the 14 N production from 
the transmutation of both 12 C and 16 via the CNO 
cycle at the bottom of the convective envelope. 

At log(£fi na i — t) = 5.08 after eleven TPs, of which the 
last six were followed by efficient dredge-up, the total 
14 N enhancement has risen to 1.83 dex and the envelope 
He abundance is Y — 0.356. While the 14 N abundance 
is now in reasonably good agreement with the estimate 
in the "extreme" abundance mixture, the value of Y is 
just reaching the lower limit of the range that has been 
inferred from the isochrones. Furthermore, the O abun- 
dance has decreased by only ~ 0.5 dex and the C abun- 
dance has increased in total by < 0.3 dex. These reduc- 
tions are both significantly less than the values that de- 
fine the "extreme" abundance distribution. However, for 
carbon, in particular, it should be kept in mind that the 
"extreme" composition is based in part on the measured 
abundances in RGB stars, w hich are known to experience 
extra -mixing processes (e.g. Denissenkov & VandenBerg 



20031 that cause the C abundance to decrease with in- 



creasing lumino sity as they climb the giant branch (|Grat- 
ton et al.||2000[ ) 



In the more advanced pulses (starting around the tenth 
TP), each dredge- up event will mix the entire intershell 
helium as well as a small layer interior to the bottom of 
the He-shell flash convection zone where X( 12 C) = 0.63 
and A( 16 0) = 0.35. The dredge-up into and below the 
former He-shell flash convection zone is shown in Fig.[5j 
The He-shell flash the convection zone combines 12 (J- 
and 16 0-rich material from deeper layers (< 0.9912 M ) 
with the 14 N-rich H-shell ashes (between 0.9924 and 
0.9927 M ). At the time that the profiles are shown the 
convection zone has not yet reached its full extent, which 
occurs when the upper boundary has approached the lo- 
cation of the former H-burning shell (0.9927 M ). At 
this point the 12 C and 16 mass-fraction abundances in 
the He-shell flash convection zone are X( 12 C) — 0.046 
and A( 16 0) = 0.012, which is considerably lower than 
the typical intershell abundances in lower mass models. 

After the end of the He-shell flash, the dredge-up (see 
the bottom panel of Fig.fsj) proceeds into the H-free core. 
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Figure 5. Profile of CNO isotopes and 4 He during the He-shell 
(approximately at the end of the first 1/3 of the duration of the He- 
shell flash convection when the convection zone is still growing out- 
ward in mass) and during the third dredge-up (when 1.3 X 10 Mq 
out of a total of eventually 1.9 X 10 — 3 Mq have been dredged-up) 
following the He-shell flash. 

The profiles illustrate how 14 N is produced in-situ as 
12 C streams through the hot convective boundary at 
m r = 0.9915 Mq. In the outward direction, it is the 
simultaneous burning of 14 N and the convective mixing 
in the envelope that causes the sloped chemical abun- 
dance profiles. The 4 He abundance just below the con- 
vective boundary is reduced from the value obtained at 
earlier times (see the top panel) due to nuclear burning, 
which forms a small additional amount of 12 C and 16 O 
to be subsequently dredged-up. However, close inspec- 
tion of the models shows that the bulk of the 12 C and 
16 O that is mixed to the surface does not come from the 
region occupied by the former convective He-shell flash, 
but originates instead from the layer just beneath. The 
deepest mass coordinate reached by the third dredge-up 
in this event is 0.9908 M . Therefore, the amount of 12 C 
coming from the region below the former Hc-shcll flash 
convection zone is 0.0003 M x 0.63 = 1.8 X 10~ 4 M Q , 
which is approximately the same as the 12 C content of 
the envelope at this point. Thus, the visible spikes in the 
12 C surface abundance evolution (Fig.[4j and much of 
the overall enhancement of CNO is due to deep dredge- 
up that reaches below the He-shell flash convection zone. 



This deep dredge-up will also bring 16 from below the 
He-shell flash convection zone into the envelope. In this 
model the dredge- up of 16 O outperforms the capacity of 
hot-bottom burning to reduce the 16 O abundance in later 
TPs. As a result, this model predicts an increasing 16 
abundance, which would eventually lead to an enhanced 

abundance in the AGB wind ejecta compared with the 
initial O abundance. This is clearly in contradiction to 
the target abundance in the "extreme" mixture. While 
this deep dredge-up will also mix all of the He in the 
helium shell into the envelope, this contribution to the 
evolution of the envelope helium abundance is negligible. 
The bulk of the He comes from hot-bottom burning dur- 
ing the interpulse phase. During the 15 fully developed 
thermal pulses (i.e., not counting the first five TPs that 
have insignificant dredge-up) the He mass-fraction abun- 
dance increases by AY = 0.02, reaching 0.372 in the last 
computed model. Therefore, a value of Y of « 0.4 could 
be obtained after another 21 thermal pulses, which is in 
all likelihood a very plausible scenario. 

In conclusion, our standard intermediate-mass AGB 
stellar model is able to generate wind ejecta with He as 
well as N abundances in the ranges required by the "ex- 
treme" mixture. As a matter of fact, the total N enhance- 
ment according to our model is 2.5 dex, which exceeds the 
enhancement specified in the "extreme" abundance dis- 
tribution. However, it is possible to alter the abundances 
of C, N and O in the models by fine tuning (reducing) 
the convective boundary mixing (CBM) parameter at the 
bottom of the convective envelope. This may be neces- 
sary, in fact, given that there will be many more thermal 
pulses in a real star than we have computed here, with 
the consequence that the envelope enrichment would be 
stretched out over a much larger number of dredge-up 
events. However, because the 12 C for the production of 

1 N during hot-bottom burning and 16 O come from the 
same region in the star, it is not possible to obtain a 
significant N enhancement and, at the same time, a sig- 
nificant O reduction (to be consistent with the "extreme" 
abundance distribution) without further modifications of 
our model assumptions. 

3.2.2. Increased mixing-length parameter 

In addition to boundary mixing uncertainties, some 
assumptions need to be made about the mixing-length 
parameter, «mlTi when modeling convection in the en- 
velopes of AGB stars. It is commonly assumed that the 
value of ajyiLTj as calibrated by fits to the solar parame- 
ters, can be applied universally to all convection zones 
during all evolutionary phases. However, both simu- 
lations as well as semi-empirical evidence suggest that 
the deep convection zones of giant envelopes, compris- 
ing essentially fully convective configurations, are better 
described by an increased mixing-length parameter. A 
variable (i.e., non-constant) «mlt parameter was already 
su ggested by the radi ation-hydrodynamics sim ulations 



( 1999 ) , also see the later work by Robin 



Bp ecifica lly relevant to giant stars, 

pool ) 



by Ludwig ct al. 
son et al'l (|2004 

Porter & Woodward (|2000j) presented 3D-simulations 
ot deep envelope convection that were best reproduced 
within the mixing-length picture if «mlt = 2.68, assum- 
ing th e formulation of the MLT given in |Cox fc Giuli| 
(19681. In addition, a semi-empirical determination of 



qmlt was obtained from the the modeling of the pulsa- 
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tion of highly-evolved, variable, inter mediate-mass stars 
(4 - 6 Mq) in the LMC and SMC by |McSaveney et al. 
( 2007 ) when they derived spectroscopic abundances, For 
three giants, they determined significantly larger values 
of Q:mt,t> ranging f rom 2.2 to 2.4 — also based on the 



m5 mlt 



Cox k, Giuli| ( |1968[ ) version of the MLT. Taken together, 
there is ample justification to explore the effects of as- 
suming a larger value of qmlt- Numerical experiments 
with enhanced values of Qmlt have been carried out for 
intermedi ate mass AGB stars a nd su per-AGB st ars pre- 



usmg 



viously by Karakas et al. (|2012|) and|Siess (|2010 
a post-processing code) 

To investigate this issue, we have computed two evo- 
lutionary sequences starting at the first thermal pulse of 
the standard sequence that was described in the previous 
section. For the convective boundary mixing (CBM) pa- 
rameters, we used f ce = 0.002 and / p dcz = 0.001, while 
at all other convective boundaries, we adopted / = 0.01. 
One sequence had «mlt = 1-73, as assumed in the stan- 
dard sequence, and the other had Omlt = 2.40. Both 
sequences were followed through the initial 5 to 6 TPs 
with no or little third dredge-up, and the computations 
were halted after another two full TP cycles with fully 
developed third dredge-ups had been completed. The 
evolutionary properties of the models over those two TP 
cycles are compared in Table|3| The «mlt = 2.40 cal- 
culations show, on average, higher interpulse H-burning, 
and peak-flash He-burning, luminosities. For the high- 
c^mlt case, the third dredge-up is more than twice as 
deep, and it reaches below the bottom of the convection 
zone boundary established by the previous He-shell flash 
at the second fully developed dredge-up, as compared 
with the standard-«MLT case wherein the dredge-up pro- 
ceeds into about two-thirds of the former He-shell convec- 
tion zone. This standard-«MLT case predicts a shallower 
third dredge-up compared with the sequence described 
in Sect. |3.2.j"1 because the adopted value of f cc is smaller. 

The convection assumptions, primarily parameterized 
through cumlt and / cc , significantly alter the evolution 
of the envelope abundances. Over the first two pulses 
with significant dredge-up, the enhancement in the He 
abundance is a factor of 6.6 higher for the cumlt = 2.40 
case than in the models for «mlt = 1.73. Also, in the 
qmlt = 1-73 case, the increase in AX( 4 He) per pulse 
at later times is consistent with the average value for 
AX( 4 He) reported for the standard case in the previous 
section. Since, in a differential sense, the same is true 
for the aMLT = 2.40 case, we find that intermediate- 
mass AGB models with a larger value of the mixing- 
length produce He from hot-bottom burning more effi- 
ciently. Insofar as the CNO elements are concerned, the 
deeper dredge-up that occurs when «mlt = 2.40 means 
that significantly more C and O is brought into the en- 
velope, even from below the former He-shell flash con- 
vection zone. However, the more efficient hot-bottom 
burning in the model with higher «mlt means that, not 
only is N produced efficiently, but even O is depleted ef- 
ficiently. Over the first two pulses of the ckmlt = 2.40 
sequence, the O abundance is depleted by ~ — 0.6dex in 
comparison with the initial abundance. 

3.2.3. Conclusion 

We conclude that, with the right combination of the 
convection model parameters q;mlt and f ce , it should be 
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Figure 6. Surface abundance evolution of 5 Mq stellar evolution 
track with enhanced mixing-length parameter ctmlt = 2.4 (top 
panel, ffl na i = 9.1641578 X 10 7 yr) and surface abundance evolution 
for 7Mq track with standard mixing assumption (bottom panel, 
tfinal =4.6642752 X 10 7 yr). 



possible, in principle, to account for He mass-fraction 
abundances up to Y ~ 0.40, as well as sufficiently 
enhanced N and depleted O, in the wind ejecta from 
intermediate-mass AGB stars. Indeed, to obtain the low 
oxygen abundances that have been observed, it seems to 
be necessary to adopt a high value of c«mlt- On the other 
hand, the increased effectiveness of third dredge-up that 
is also found when a larger value of «mlt is assumed, 
may lead to an enhancement of N that is too large for a 
given increase in Y. We have therefore constructed one 
additional 5 M stellar model sequence with «mlt = 2.4, 
but with a further reduction in the value of the / ce pa- 
rameter. This case has been evolved over ten thermal 
pulse cycles with fully developed dredge-up (see the top 
panel of Fig. [6]). While the N enhancement now reaches 
~ 2 dex and the He abundance reaches the average value 
of the range that is indicated by fitting isochrones to 
the bMS, the O abundance has been depleted by about 
1 dex, just as required by the "extreme" abundance dis- 
tribution. However, even this model sequence does not 
show a significant C depletion. 

A real intermediate-mass AGB star will likely do more 
thermal pulses than we were able to compute here. This 
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means either that the real star may have convection prop- 
erties that will stretch out the chemical abundance mod- 
ifications seen in our stellar models over more thermal 
pulses, or that the ejecta of the AGB star are more ex- 
treme in their abundance patterns than indicated by the 
"extreme" metals mixture, which would leave some room 
for the dilution of the ejected winds with unprocessed 
"standard" abundances before forming the bMS genera- 
tion of stars. As discussed earlier the "extreme" abun- 
dance mix may in fact be an upper limit for the bMS in 
ui Cen, and a somewhat less extreme abundance mix in 
this cluster would add to the possibility of dilution. 

With these caveats in mind, it seems very reasonable 
to conclude that intermediate-mass AGB star models are 
capable of producing abundance patterns very similar 
to that listed in Table[l] under the column heading "ex- 
treme" . 
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3.3. Super- AGB model 

For somewhat higher initial masses, stars have a core 
mass at the first thermal pulse which is larger than ap- 
proximately 1.05 Mq, in which case, C bur ning is ignited 
and those sta rs will form a ONeMg core ( Garcia-Berro 
fc Iben|[T994] ). These are the so-called l super-AGB" 
stars, which, like their intermediate-mass AGB star coun- 
terparts with CO cores, have periodic thermal pulses, 
dredge-up events, and hot-bottom burning. It has been 
suggested that these super- AGB stars could provide the 
chemically peculiar ejecta that are required in order for 
a stellar population like the bMS in u> Cen to form. 

3.3.1. Deep second dredge-up mixing in the 7Mq case 

We have calculated stellar evolution sequences with an 
initial mass of 7 M that adopt the same initial abun- 
dances, and employ the same physical assumptions, as in 
the 5 M sequence discussed above, but explore uncer- 
tainties in convection assumptions. 

In Fig.|6] we show the evolution of the surface abun- 
dances for our standard c ase. As reported for the 
lower mass models in Sect. |3.2TT| the chemical compo- 
sition changes as a result of the first dredge-up (at 
log(tfi na i — t) — 6.7), while the second dredge-up begins 
to have an impact at log(tfi na i — t) = 5.1. The He abun- 
dance gradually increases as the bottom of the convec- 
tive envelope penetrates into the ashes of the H-burning 
shell. In some super- AGB models, depending on mass 
(see Sect. 3.3.2 1 and details of the physics assumptions, a 
He-shell hash convection zone has developed before the 
end of the second dredge-up. The quenching of the flash 
causes a fur ther dredge-up, also refered to as a "dredge 



out" event (Iben et al. 1997 Ritossa et al. 1999 Siess 



20071. This final mixing episode is shown in more detail 



m Fig.[7j Depending on the mass the second dredge-up 
proceeds into the He-free core. This mixing into the He- 
free core adds additional C, and later on also O to the en- 
velope, as seen in Fig.[7|and Fig.|6]at log(ifi na i — t) = 4.25. 
This dredge-out mixing increases the envelope O abun- 
dance in our 7 M model by ~ 0.77 dex, which is com- 
pletely at odds with the ~ 1 dex reduction that is re- 
quired in order to match the abundance given in the 
target "extreme" mixture. The C abundance is also sig- 
nificantly enhanced, by 1.6 dex, with respect to the intial 
value. 
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Figure 7. Top panel: detail of the last phase of the second dredge- 
up of the Mzams=7M0 case; bottom panel: simultaneous surdace 
abundance evolution. After a so-called "dredge-out" phase the 
bottom of the core convection zone continues to penetrate (on La- 
grangian coordinate) into the He-free core, bringing first 12 C and 
then le O into the envelope. The quantitative details of these mix- 
ing processes are sensitve to uncertain physics assumptions (see 
text). H- and He- free boundaries correspond to the mass coordi- 
nate where the mass fraction is respectively < 10 — 4 . 

One may think that efficient hot-bottom burning in 
super-AGB stars, possibly boosted by the adoption of a 
higher value of «mlt could reduce even this additionally 
dredged-out O and C. Although individual third dredge- 
up events cannot be identified in the plot which shows 
the evolution of the envelope abundances, the 7M Q se- 
quence has been evolved through 18 thermal pulses with 
deep dredge-up. They cause the helium abundance to in- 
crease by AA( 4 He) = 0.0002 per TP, which is the same 
as that found in the standard MLT sequence presented 
in Table[3] In order to reach a He mass-fraction abun- 
dance of 0.40 another 105 thermal pulse cycles would be 
needed, which is a very likely possibility considering the 
short interpulse periods and a reasonable range of possi- 
ble mass-loss rates. 

During the same initial 18 thermal pulses, the O abun- 
dance decreases by 0.16 dex; consequently, over another 
~ 100 TPs, a significant O depletion could be achieved. 
However, even if hot-bottom burning is able to further 
increase the abundance of He, and reduce the O abun- 
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Figure 8. Lifetimes of intermediate-mass AGB and super-AGB 
for the "standard" initial composition mix given in T able]!] deter- 
mined from a grid of stellar models described in Sect, |3,3^| 

dance, by the desired amounts, the N abundance will 
become far too large. At the end of 18 thermal pulses, 
the 14 N abundance has already increased by more than 
2.3 dex and futher depletions of O will only add to the 
production of N. While the total CO enrichment in CO- 
core AGB stars can be controlled by fine-tuning the con- 
vection parameters, this is not obviously possible given 
the substantial enrichment of some super-AGB envelopes 
during the deep second dredge-up that reaches into the 
He-free core. 

3.3.2. Initial mass dependence 

The deep second dredge-up depends on assumptions 
of convective boundary mixing a well as on stellar mass. 
In order to determine the mass dependence of this mix- 
ing, as well as the upper and lower limiting masses for 
the super-AGB phases and the lifetimes for intermediate- 
mass AGB and super-AGB stars, we have calculated a 
fine grid of models with 4.0 < Mzams/Mq < 8.80 with 
AMzams — 0.2 M Q . Initial abundances are again ac- 
cording to the standard mix, and convective boundary 
mixing for H- and He-core burning with / = 0.014 was 
included. At all other convective boundaries / = 0.010 
was assumed. For this grid v. 3372 of the MESA st ellar 
evolution code w as used, together with the NACRE ( An- 
gulo et al~||1999 ) reaction rates. From these calculations 
we find the lifetimes shown in Fig. [8] 

For Mzams > 5.6 M Q with H-free core masses of Mh > 
1.05 M Q the models show off-center C ignition, which in 
these models does not burn to completion in the central 
region in most cases. Deep second dredge-up that pro- 
ceeds into the He-free core is found for M Z ams > 6.8 M , 
and in these models the maximum He-free core mass 
before the reduction through the second dredge-up is 
-^He > 1.245 M Q . Above this mass models predict large 
dredge- up of primary CNO as described in Sect. |3. 3.1~] 
For M Z ams > 7.6 — 7.8 M Q the core mass after the sec- 
ond dredge- up is at or just below the Chandraskhar mass, 
and more massive stars are expected to explode as super- 



3.3.3. Convective boundary mixing effects 



Table 3 

Comparison of two thermal pulse cycles with standard 
(c*MLT = 1-73) and increased mixing-length parameter 
( a MLT = 2.40): difference in H-shcll intcrpulse and He-flash peak 

luminosities, dredge-up depth in Mq for second TP with deep 
dredge-up discussed in the text, surface He mass fraction increase 
over the two TP cycles as well as change of surface C, N and O 
over the two considered TP cycles. 





£*MLT = l-73 


«MLT=2.40 


A log L H [2.40 - 1.73] 




0.2 


AlogL Hc [2.40- 1.73] 




0.3 - 0.5 


Am dup / Mq 


7.1 x 10" 4 


1.66 x 10~ 3 


AX( 4 He) 


0.0004 


0.0023 


A log X (C) 


0.28 


0.50 


AlogX(Af) 


0.26 


0.50 


AlogX(O) 


-0.10 


-0.51 



We have calculated another 7M sequence without 
convective boundary mixing at the bottom of the con- 
vective envelope, but its behaviour with regard to the 
yields of He and CNO over the first dozen thermal pulses 
is very similar to the standard case. It shows the same 
deep second dredge- up. F or these more massive super- 



AGB stars (cf. Sect. 3.3.2) the yields are dominated by 
the second dredge- up. The third dredge-up will only dou- 
ble the C+N+O abundance over the 18 computed ther- 
mal pulses (in the form of N) , while the sequence without 
convective boundary mixing will not add additional CNO 
because there is no third dredge-up. 

Indeed, significant differences exist when convective 
boundary mixing is taken into account or when it is ne- 
glected altogether, in that the former is accompanied by 
effective third dredge- up with A < 1, while in the lat- 
ter instance, the third dredge-up is absent. This has the 
important consequences that the cores of models with 
CBM do not grow, and they are therefore incapable of 
reaching the Chandrasekhar mass. As a result, they will 
become white dwarfs instead of evolving into ONeMg 
core-collapse supernova — unless they reside in a binary 
system, in which case they may still explode through an 



accretion-induced collapse (Fryer et al. 1999). 



Another, potentially even more relevant consequence is 
the presence of z-process-condition s in super-AGB mod- 
els with CBM ( |Herwig et aT]|2012[ ). Briefly, we find that 
dredge-up may - due to the short thermal time-scale at 
the core-envelope boundary - enter into the intershell 
immediately after the peak He-burning luminosity has 
been reached, and the He-intershell region is still con- 
vectively unstable. Then, depending on the convective 
mixing assumptions, protons may be convectively mixed 
into He-shell flash conditions, a situation in which high 
neutron densities of N n ~ 10 15 cm~ 3 would be present 



dCowan fc Rose|1977| [Campbell et al.|2010||Herwig et al" 

201 IB. The details of the i process in super-AGB start 



process m super-AGri stars 
are still very uncertain, mostly due to our inability do 
correctly model the interaction between convective mix- 
ing and rapid, combustion- like burning of 12 C + protons 
with the spherically symmetric assumption of stellar evo- 
lution , and 3D simulations, such as those in Herwig et al. 
(2011 1 wi ll be required. Howev er, within the present un- 



certainty Herwig et al. (20121 do consider a significant 
production of many n-capture elements in super-AGB 
stars from i process, including La, a possibility. The im- 
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plicat ions of this possibilty will be discussed again later 
(Sect.[K2|. 



3.3.4. Conclusion 

Super-AGB stars with lower initial masses than ~ 
6.8 M Q and consequently lower He-free core masses (< 
1.245 Mq), would not suffer from the additional C and O 
envelope enrichment from deep second dredge-up. Stars 
in this mass range would in principle behave sim ilar to 



the massive (CO core) AGB stars described in Sect. 3.2.1 
and could have ejecta with the abundance pattern of the 
extreme distribution. The contribution of the most mas- 
sive super- AGB stars must be minor due to the unavoid- 
able and excessive O and N abundances that would char- 
acterize their ejecta. In order to identify the mass ranges 
and exact outcomes quantitatively in a more reliable way, 
the properties of convection in the deep stellar interior 
must be investigated in more detail. 

In addition to the uncertainties in the treatment of 
convection, mass loss assumptions will significantly effect 
final and quantitative abun dance predictions for ejecta 
( Ventura &: D'Antona|20lT |. In our models we have kept 
mass loss low in order to study the impact of convection 
uncertainties in isolation. 

4. SUMMARY AND DISCUSSION 

In this section we provide a brief summary of our re- 
sults, and then suggest a scenario for the formation of 
discrete populations of stars in uj Cen, with an emphasis 
on the bMS. 

4.1. The main results of this investigation 

Consistent with previous findings, we have found that 
the location of the bMS in the CMD of uj Cen can be 
reproduced by stellar models only if they assume sig- 
nificantly enhanced He abundances, in the range Y = 
0.35 — 0.40. Spectroscopic observations for giants in 
this system, and for one blue main-sequence star in 
NGC2808, suggest that it is not only helium but also 
the CNO elements that have peculiar abundances, with 
mean values of [m/Fe] that have been described as the 
"extreme" Pop. II metals mixture (see Tablcll}. The 
bMS stars appear to have a higher iron content hy about 
a factor of two than the dominant more metal-poor stel- 
lar population, which has [Fe/H] rj —1.7. 

The other, equally puzzling, CMD feature is the red- 
dest MS that is commonly referred to as MS-a. We 
have shown that it is possible to match the photome- 
try of these stars remarkably well using isochrones for 
the measured [Fe/H] value (ss —0.7), provided that they 
also have rather high helium and CNO abundances. In- 
deed, it is possible that the bMS and MS-a populations 
differ only in terms of their iron abundances; i.e., that 
they are characterized by the same high values of Y and 
the "extreme" [m/Fe] ratios. If this suggestion is correct, 
then all of the stellar populations in uj Cen would appear 
to have similar ages (within ss 1 Gyr given the chemical 
abundance uncertainties). For the MS-a stars to be sig- 
nificantly younger than, in particular, the most metal- 
deficient stars, the former would need to have much 
higher [CNO/Fe] abundances than those implied by the 
"extreme" mixture, since the latter abundances coupled 
with an old age (13 Gyr) are already needed to explain 
the very faint subgiant branch of the MS-a component. 



We have further established, through simulations of 
AGB and super-AGB stars, that both of these types of 
stars, when forming out of the abundance mix which has 
been assumed for the first and dominant [Fe/H] = — 1.7 
population, are capable of ejecting winds with the high 
helium abundances that are needed to explain the CMD 
locations of the bMS stars. However, for the source of 
this material, we favour the massive AGB stars with CO 
cores and the less massive super-AGB stars with ONeMg 
cores. The wind ejecta of the most massive super-AGB 
stars with core masses > 1.245 M Q are likely too enriched 
with N relative to He. In any case, even though either 
of the AGB models can eject matter with sufficiently 
high He abundances, it is not clear that they can pro- 
duce ejecta with much larger He abundances than the 
upper limit that has been set from the CMD of the bMS 
(Y < 0.4). This imposes limits on the acceptable amount 
of dilution of the AGB material with gas that may still 
be present in the cluster or accreted from the interclus- 
ter medium in order to allow bMS stars to form with the 
indicated high He abundance pattern. Similarly, stellar 
populations with very low O abundances set t ight con- 
straints on the amount of allowable dilution (D'Ercole 



et al. 2011a 



The obvious question is then how to isolate the AGB 
ejecta from other gas that could be present; i.e., how 
to form the bMS stars out of just the ejecta from mas- 
sive AGB stars instead of, for instance, the C-rich ejecta 
of somewhat lower-mass AGB stars (or possibly the Fe- 
rich ejecta of the more massive stars that explode as su- 
pernovae). The observationally established central con- 
centration of stars with the "extreme" abundance mix 



(Johnson & Pilachowski 2010) may very well be an im- 
portant clue. 

4.2. A scenario for the formation of multiple 
populations in globular clusters 

We need to combine the stellar evolution results de- 
scribed thusfar with additional information in order to 
propose a general scenario for the formation of multiple 
populations, such as the bMS and MS-a/RG-a compo- 
nents of uj Cen. In particular, it is necessary to consider 
the early evolution of globular cluster progenitors in the 
gravitational potential of the host galaxy in combination 
with the cooling properties of AGB ejecta. To be sure, 
several aspects of this scenario have already been sug- 
gested and elaborated upon elsewhere, as pointed out 
below. 

4.2.1. The progenitor of uj Cen and the orbit in the Galactic 

potential 

The present-day orbit of uj Cen consists of a super- 
position of several oscillation components (Fig. [9]), which 
suggests a complex capture mechanism that may very 
well have involved additional merging components. Ac- 
cording to the backwards in time integrated orbit uj Cen 
has passed through the Galactic disk 13 times over the 
past 5 x 10 s yr (Fig.[9| , implying an average time in- 
terval between Galactic plane passage events of ~ 40 
million years, with some variance. This estimate is ap- 
proximate because of the uncertainty associated with the 
present-day velocity of uj Cen, and furthermore, the time 
interval may have shortened significantly since the cap- 
ture of uj Cen's progenitor, which may have been a dwarf 
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Time/Gyr [0 = presenL] 

Figure 9. The orbit of LoCen integ rated assuming the pre sent-day 
coordinates and positions given b y Dincscu ct al. ( 19991 and the 
Galaxy model described inlMichel-Dansac et al. 1120111. The orbit 
is computed by integrating backwards in time (i.e., alter changing 
the sign of the present-day velocity) and shown for the past 500 
Myrs. The solar circle is indicated, at Rq = 8.5 kpc, by a dotted 
line. Note that uCen crosses the disk of the Galaxy roughly every 
30-50 Myrs. 

galaxy (iBekki fe Freeman|2003[ |B5ker|2008| iMarcolini fe 
DTjcole||2008| ). 

In tins primordia l scenario ol globular cluster forma- 
tion ( |Padoan et al. [1997 ), w Cen would have been much 
more massive (e.g ., by a factor of ~ 25 according to Bekki 
& Freemanl 2003 ) when it was accreted, and it would 



have lost substantial fractions of its stellar and halo mass 
through tidal shocks at each pericentric passage. This 
process has been demonst rated through, for examp le, 
the n-body simulations by |Peharrubia et aL ( 2008 ) of 
dSph galaxies orbiting in tne potential ol galaxies like 
the Milky Way. These computations show that, depend- 
ing on the orbital parameters (specifically the apocen- 
tric/pericentric ratio), within the first couple of pericen- 
tric passages (i) ~ 90% of the dark matter halo is lost 
before any significant loss of stars occurs and (ii) mass is 
lost through outside-in "onion-peel" stripping. Through 
such tidal mass loss and evaporation episodes, it is quite 
conceivable that a much more massive progenitor would 
have ended up in uj Cen's present state. 

For globular clusters, Galactic plane passages have long 
been associated with the ram p ressure stripping of th e 
gas content of globular clusters (Tayler & Wood 1975). 
We therefore assume that dark matter, stellar and gas 
components of the merging progenitor of uj Cen will be 
lost in varying proportions (as have , e.g. |Valcarce & 



Catelan||2011| |D'Antona et al.||2011b[ ), that will depend 
on the exact mass and orbit as well as the state of the 
Galaxy at the time of each interaction. The outcome 
would in any case be a complete purging of the gas in 
Galactic plane passages. 

4.2.2. Supernova feedback and purging 

Studies of supernova feedback in the context of the 
formation and evolution of galaxies show that star for- 
mation correlates with the mass of the galax y and anti- 
correlates with the energy of the supernova ( |Scannapieco 



et al. 2008). Supernova feedback in low-mass systems 
leads to more severely limited and bursty star forma- 
tion. If the progenitor of uj Cen had low enough mass in 
its early history, supernova feedback could have caused 
the purging of all of its gas. Indeed, an alternative to the 
Galactic plane passage purging is gas purging through su- 
pernova from sta rs with initial masses that bracket the 



AGB mass range (D'Ercole et al. 2008 Conroy & Spergel 
201 1 1 : stars more massive than super- AGB stars explode 



as type II supernovae, while those near the lower mass 
limit for massive AGB stars may lead to prompt (single- 
degenerate) supernovae of type la. In each of these cases, 
depending on the mass of the stellar system, SN explo- 
sions could purge the gas from it. 

The cluster would, however, somehow need to be able 
to retain some of the iron from SN (la ?) explosions 
at later times in its evolution in order to enrich the gas 
out of which the relatively metal-rich MS-a population 



forme d (iPancino et al.| [2002; Ori glia et al.|2003 Pancino 
et al. 20ll| ). Similarly, at the epoch when stars were 



evolving into massive AGB stars, the potential well of 
uj Cen likely has to be deep enough for this system to 
retain the Fe ejecta from the small range of initial masses 
corresponding to the lowest mass type II SN, which have 
the lowest explosion energies. Most of the Fe that was 
produced by SN from more massive progenitors would 
have been blown out of the GC along with any intra- 
cl uster gasj^] 

D'Ercole et al. (2008) have pointed out that, if su- 
pernova from the electron-capture core collapses of the 
ONeMg cores of sup er- AGB stars have si gnificantly lower 
explosion energies (Dessart et al. 2006), the (possibly) 
Fe-rich ejecta of those SIN may be more easily retained 
by the cluster. This may be a viable alternative expla- 
nation of the small increase of the iron abundance in the 
bMS stars compared to the measured [Fe /H] of sa —1.7 in 
the do minant initial stellar population. Marcolini et al.| 
(2007) have investigated such aspects of the evolution of 
uj Cen through 3D hydrodynamic simulations of an iso- 
lated progenitor system that was assumed to be a dSph 
galaxy. They emphasized the differences between the SN 
II and la effects on the gas polution and the possibility of 
inhomogeneous enrichment of gas with SN la ejecta, but 
they admitted that tidal interactions with a host galaxy 
may have additional important consequences. 

In any case, there are two, possibly complementary, 
processes available for the purging of all gas from the 
cluster — SN purging or periodic ram pressure/tidal 
shock purging through Galactic plane or pericentric pas- 
sages. Which of these mechanisms dominates will de- 
pend o n the d etailed dynamical evolution of the system. 
In Sect. 4.2.4 we will explain why we favor the purging 
of gas via Galactic plane passages in the case of uj Cen 
because of the unique effects that this process can poten- 
tially have on the star formation and chemical evolution 
histories of such a system. 

9 In this picture, the cluster would have also expelled most of 
the wind materi al from rapidly ro tating massive stars, which was 
proposed by |Decressin et al.| l |2007| l to be the production site of the 
He-enriched material that went into the formation of such extreme 
po pulations as the bMS in 10 Cen. However, it has been shown 
by |Romano et ah] l |2007^ that these special massive stars cannot, 
by themselves, account tor the required abundance patterns of the 
bMS. 
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Figure 10. The radiative cooling coefficient for a gas having the 
indicated chemical abundances (also see the text). The loci with 
filled circles superimposed represent the sum of the contributions 
to A/(n e nH) from H, He, C, N, O, and Ne, as well as free-free 
emission (which dominates at logT <^ 5.8). For a sake of clarity, 
only a few of the contributions to the total cooling curves from 
individual elements have been plotted. 

4.2.3. The Cooling of the Gas Ejected by AGB Stars 

Cooling flows from AG B ejecta have been investigated 



by D'Ercole et al. (2008), on the assum ption of a cool- 



ing function tor the solar metallicity by Rosen & Breg- 



man 



(1995). However, the ejecta of the first generation 
AGB stars have a much lower metallicity, and therefore 
one would expect that the cooling of this material would 
be much less efficient. Since the cooling flow of AGB 
ejecta to the center of a system like uj Cen is an impor- 
tant ingredient in most scenarios (ours, and those by, 
e.g.JD'Ercole et al.||2008l ID'Antona et~aLl|2011a[ [BekEl 



2011 1 for t he origin of the He-rich second generation (see 



4.2.4 ), it is worthwhile to examine the cooling prop- 



Sect _ 

erties of a gas having the "extreme" metals mixture, and 
how it differs in comparison with that predicted for for 
"standard" (i.e., typical) Pop. II abundances. 

While the actual abundances in the uj Cen bMS stars 
may be somewhat milder than in our "extreme" abun- 
dance mix, those abundances out of which the stars 
formed may have b een subject to some limited amount 
of dilution (cf. Sect. 4.1 ). Any assumption of dilution re- 
quires the AGB ejecta to be more extreme than the bMS 
stars in uj Cen. What exactly the abundance of the cool- 
ing flow is will depend on when the dilution takes place, 
something that is difficult to know. In any case, it is ap- 
propriate to consider the cooling properties of material 
with "extreme" composition as representative of what the 
AGB ejecta may look like. 

In the limit where the photoionization caused by the 
stellar radiation field can be neglected; i.e., when col- 
lisional ionization conditions apply, the radiative cool- 
ing coefficient (in ergcm _3 s _1 ) va ries as a function of 
log T approximately as shown in Fig[l0]for three different 
heavy-element mixtures. This plot was generated using 
the computer code developed by VandenBerg (1978) to 
model the outflows of gas from present-day globular clus- 
ters to try to explain why the gas that is ejected into the 
interstellar medium (ISM) through normal (low-velocity) 



mass-loss processes is not observed. (When these sys- 
tems pass through the Galactic disk, any gas that might 
have collected is expected to be swept out by the ram 
pressure exerted by the local ISM, but gas should accu- 
mulate to detectable levels over the rest of their orbits 
in at least the most massive systems.) If photoionization 
is unimportant, then the ionization state of the gas can 
be determined simply from the condition that the num- 
ber of collisional ionizations per unit time for any ion is 
exactly equal to rate of electron recombinations, which 
is independent of the density. Once the ionization state 
has been calculated, the radiative energy losses due to 
free-free and free-bound transitions, and from the tran- 
sitions to lower energy levels from high-excitation states 
that were populated by recombination or inelastic elec- 
tron collisions, can be determined (see the VandenBerg 
paper for a detailed description of the relevant physics). 

Since the energy loss rates due to recombination radia- 
tion and collision-induced line emission are proportional 
to n e nz,z, where n e is the electron density and nz, z is 
the number density of atoms having an atomic number 
Z and an ionic charge z, the radiative cooling coefficient 
A/(n e nu) can be derived for any mixture once the frac- 
tional ionization, nz tZ /nz, has been calculated and the 
abundances of each element are specified relative to hy- 
drogen (i.e., Uz/tih). The solid curves in Fig. 10 as- 



sume the solar metal abundances given by |Asplund et ak] 
( |2009| , while the dot-dashed and dashed curves give, in 
turn, the temperatu re dependence of A/(n e nu) for an 
Asplund et al. (2009) mixture having [m/H] = 0.4 for all 



of the a-elements ( standard" composition in TabletO) 
and then scaled to [Fe/H] = —1.4 and for the "extreme^ 
mix from Tablejl] scaled to [Fe/H] = —1.4. In these latter 
cases, the adopted mass fractions of helium correspond 
to Y = 0.25 and Y — 0.38, respectively. (Since these 
results may be of some use to future investigations that 
take into account the cooling of the gas, the total cooling 
coefficients, as represented by the loci with filled circles 
attached to them, have been listed as a function of log T 
in Table g) 

For the sake of clarity, just a few of the contributions to 
the total cooling rates, due to the individual metals, are 
shown. Nitrogen is of particular interest in view of the 
wide range in its abundance in the three mixtures, and 
it is quite evident that the cooling due to this element 
in the "extreme" [Fe/H] = — 1.4 mixture is much greater 
than that predicted for the other two cases (including, in 
particular, solar abundances). It is also much greater, by 
factors of ^ 60, than the contributions due to C and O 
because the latter are much less abundant. As helium is 
mainly responsible for the bump in the radiative cooling 
coefficient at log T ~ 4.9, the effects of increasing Y 
from 0.25, in the "normal" mix, to 0.38, in the "extre me" 
mix, are clearly significant. The main point of Fig. [10] 
is that a gas with high N and He abundances will cool 
much more efficiently than a gas having a more normal 
mix of metals and a helium abundance that is closer to 
the primordial value. Indeed, the total radiative cooling 
coefficient for the former approaches that for a gas having 
solar abundances (at least at logT J5 5.3). 

The presence of an intense stellar radiation field, such 
as that produced by extremely hot (~ 50, 000 K) uv- 
bright stars, would heat the gas and ionize many of the 
atoms and ions that are responsible for radiative cooling, 
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Table 4 

Radiative Cooling Coefficient: log(A/n e njj) 





[rc/ rlj — U.U 


r lv, / u 1 
[te/HJ 


1 A 

— — l.^t 


log J 




Standard Mix 


H/Xtreme 1V1 


3.80 


-24.694 


-24.694 


-24.694 


3.85 


-24.626 


-24.626 


-24.626 


3.90 


-24.425 


-24.425 


-24.425 


3.95 


-24.012 


-24.013 


-24.013 


4.00 


-23.489 


-23.490 


-23.490 


4.05 


-22.972 


-22.973 


-22.973 


4.10 


-22.513 


-22.514 


-22.514 


4.15 


-22.162 


-22.163 


-22.163 


4.20 


-21.992 


-21.994 


-21.993 


4.25 


-22.005 


-22.009 


-22.009 


4.30 


-22.110 


-22.121 


-22.118 


4.35 


-22.228 


-22.257 


-22.248 


4.40 


-22.319 


-22.386 


-22.360 


4.45 


-22.359 


-22.493 


-22.433 


4.50 


-22.353 


-22.589 


-22.482 


4.55 


-22.308 


-22.679 


-22.514 


4.60 


-22.232 


-22.752 


-22.509 


4.65 


-22.134 


-22.783 


-22.453 


4.70 


-22.022 


-22.736 


-22.337 


4.75 


-21.896 


-22.589 


-22.168 


4.80 


-21.761 


-22.377 


-21.968 


4.85 


-21.631 


-22.183 


-21.788 


4.90 


-21.530 


-22.079 


-21.682 


4.95 


-21.466 


-22.076 


-21.651 


5.00 


-21.443 


-22.130 


-21.653 


5.05 


-21.469 


-22.202 


-21.654 


5.10 


-21.493 


-22.266 


-21.644 


5.15 


-21.482 


-22.308 


-21.629 


5.20 


-21.454 


-22.328 


-21.632 


5.25 


-21.429 


-22.338 


-21.713 


5.30 


-21.415 


-22.342 


-21.912 


5.35 


-21.402 


-22.343 


-22.161 


5.40 


-21.412 


-22.362 


-22.381 


5.45 


-21.496 


-22.443 


-22.559 


5.50 


-21.680 


-22.602 


-22.712 


5.55 


-21.889 


-22.771 


-22.836 


5.60 


-22.046 


-22.892 


-22.931 


5.65 


-22.140 


-22.968 


-23.002 


5.70 


-22.205 


-23.022 


-23.060 


5.75 


-22.298 


-23.089 


-23.120 


5.80 


-22.467 


-23.190 


-23.187 


5.85 


-22.677 


-23.293 


-23.244 


5.90 


-22.857 


-23.363 


-23.275 


5.95 


-22.974 


-23.402 


-23.280 


6.00 


-23.031 


-23.419 


-23.266 


6.05 


-23.050 


-23.423 


-23.242 


6.10 


-23.046 


-23.417 


-23.215 


6.15 


-23.028 


-23.405 


-23.190 


6.20 


-23.001 


-23.388 


-23.175 


6.25 


-22.972 


-23.368 


-23.173 


6.30 


-22.947 


-23.348 


-23.183 



and thereby hinder or prevent the infall of gas towards 
the cluster center. However, such stars, which represent a 
much more significant source of heating than large pop- 
ulations of blue HB stars (see VandenBerg 1978), are 
quite rare (;$ 1 star is expected in present-day GCs at 
any given time although this number is expected to have 
been larger in the past). Consequently, it is difficult to 
describe the interplay between the radiation field and the 
material shed by AGB stars, without performing detailed 
simulations of globular clusters when they were younger 
and much more massive than they are now. 

It will be interesting to see what impact these more 
realistic cooling properties will have on simulations that 
include the cooling flo w of AGB stars (as we ll as that of 
primordial matter, e.g. D'Ercole et al. 20081. Consider- 
ing the "extreme" abundance mix as an upper limit of 



the boost in cooling efficiency over the "standard" mix, 
the cooling efficiency of realistic AGB ejecta are smaller 
by a factor > 4 compared to solar cooling curves employd 
be those simulations. In any case, we do assume in the 
following that AGB eject a will cool a nd accumulate in 
the cluster center (see also Bekki|20lT who comes to the 
same conclusion). 

4.2.4. A possible scenario 

The lifetime of the 7 and 5M Q stellar evolutionary 
models is 4.6 and 9.1 x 10 7 yr, respectively (Fig. [8]). The 
similarity of the time interval between successive Galac- 
tic plane passages of u Cen and the time it would take 
super-AGB and massive AGB stars to expel their He- 
and N-enriched and O-depleted wind ejecta could very 
well conspire in such a way that, by the time that the 
massive- and super-AGB stars are about to eject their en- 
velopes, the wind material from more massive stars has 
been entirely, or mostly, cleared out of the system and 
the AGB winds are released into an empty, or nearly 
empty, cluster. Stars with an extreme abundance mix 
would subsequently form, but the star formation process 
would be terminated when the remaining gas is cleared 
out during the next passage through the plane. This 
would isolate the AGB yields and enable formation out 
of pure AGB ejecta. 

The scenario for the origin of multiple populations 
starts with the progenitor evolution of u> Cen . The proto- 
cluster, possibly a dwarf galaxy (Sect. 4.2.1 1, would have 
a fully populated initial-mass function (IMF) when it 
merges with the Galaxy. Star formation may have going 
on at a low level since its formation, corresponding to 
the low-metallicity tail below the dominant (first) gen- 
eration at [Fe/H]~ —1.7 in u Cen. But somehow, we 
assume that the dominant stellar content of the ui Cen 
progenitor would be the first generation, and it would 
have been produced by a star formation burst, maybe 
triggered by the tidal interactions associated with the 
complex capture proces ses indi cated by the present-day 
orbit (Fig.|9j Sect . 14. 2. 1| see also |Bekki fc Freeman|2003| . 
Assuming that u> Cen is then in an orbit fairly similar to 
the one shown in Fig.[9j the first Galactic plane passage 
gas purging could occur ;$ 40 million years later. This 
would remove all of the remnant gas ejecta from stars 
more massive than super-AGB stars. 

Over the next 50-90 million years, until the next 
Galactic plane passage, AGB stars and possibly super- 
AGB stars would eject winds into the empty cluster, with 
abundance patterns that have helium and oxygen abun- 
dances that closely resemble those attributed to the bMS 
stars in w Cen (Sect. [3]). These ejecta would initially have 
the same spatial distribution througout the cluster as the 
AGB donor stars, but radiative cooling processes would 
cause the gas to sink preferentiall y tow ards the center of 
the cluster, as d isccuss e d in S ect. 4.2.3| 

According to Bekki (2011), the accumulation of gas 
and the subsequent formation of the second generation of 
stars will be fully under way after ~ 3 million years, and 
be largely completed in the central core within ^ 13 mil- 
lion years. This leaves enough time until the next Galac- 
tic plane passage to fully populate a second-generation 
IMF. 

Bekki suggests that the second-generation IMF may 
be bottom heavy in order to explain the high ratio of 
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second- to fir st-generation stars. However, as mentioned 
in Sect. 4.2.1 and noted by Bekki, among others, further 
mass is lost from the system through subsequent peri- 
centric passag es from the outside-in (see the discussion 
.4.2.1). In fact, the fraction of first generation 



Sect. 



stars that lias to be lost (from the outside-in) in order 
to match the observed bMS star fraction can be easily 
estimated for a given mass range of intermediate mass 
stars that release their winds between Galactic plane pas- 
sages. This mass range will of course depend on the (not 
accurately known) orbital parameters of u Cen at that 
time. Based on the calculated model grid described in 
Sect. |3.3T2] we adopt a lower limit for this mass range as 
M m i n = 4-5 M© by the requirement to produce large He 
and N enhancements, and M max = 6.8 - 7.7 M taking 
into account possibly prohibitive dredge-out (Sect. 3.3.1 1 
and excluding the high-mass regime of core-collapse su- 
pe rnova. With t his mass range assumptions the IMF 
of Kroupa (2001) implies that intermediate mass stars 
will eject 5.2 ± 2.0% of the mass of the first generation 
stars (assuming 20% of the initial stellar mass will end 
up in the white dwarf). The number of bMS stars has 
been projected to be 24 - 35% of the first ge neration 



stars flBedin et al. 2004 



Sollima et al. 2007b). These 



numbers imply that m order to account for the observed 
fraction of bMS stars 78.3 ± 8.3 - 85.1 ± 5.7% of first 
generation stars must have been lost from the outside 
(assuming a normal IM F). Th is seems to be reasonable 
in our scenario (cf. Sect. 4.2.1 ), considering also that the 
central location of the bMS population implies that it 
contributes only marginally to the mass that is lost from 
the cluster through dynamical interactions. We therefore 
assume that the second generation IMF is normal. 

Then, the second generation will also include massive 
stars. The associated supernova will explode in the very 
center of the still relatively massive cluster, which is the 
most favorable locatio n for the retention of SN ejecta 
(Marcolini et al. 20071. Depending on when these SN 



occur with, respect to the timing of the second Galactic 
plane passage, the ejecta from the most massive second- 
generation supernova may be lost during the passage 
through the disk, but afterwards, it seems plausible that 
Fe-rich ejecta from the lowest-mass SNe could be retained 
in the core region. Further down the IMF, super-AGB 
stars and possibly the most massive AGB stars with 
CO core could again contribute to what would be the 
third generation. This would again form in the cluster 
core before the third Galactic plane passage would ter- 
minate this last major star formation episode. This third 
generation would have a substantially higher iron abun- 
dance from second-generation SN, compared to that of 
the bMS, and it would have enhanced He and CNO. Fur- 
thermore, this third generation, which may be the RGB- 
a/MS-a population discussed in Sect.[2j would consist of 
significantly fewer stars than the second generation. As 
before, this star formation event would end when the 
cluster makes its fourth passage through the Galactic 
plane and the remaining gas is purged. After that di- 
minishing numbers would make it impossible to discern 
further populations in the cluster core. 

Star formation would not be impossible in other parts 
of the cluster after this fourth Galactic plane passage, but 
it is possible that, by this time, the mass of the cluster 
has decreased to the extent that it is no longer possible 



to retain the ejecta from supernovae. In this case, the 
total cluster wind ejecta may now, > 2 x 10 8 yr after 
the initial star formation burst, be dominated by the in- 
creasing number of intermediate-mass (~ 3 — 4M Q ) first- 
generation stars that are ending their lives. The ejecta 
would again be s-process rich, but also C-rich. Stars 
forming out of these "late" wind ejecta would have rather 
low [Fe/H] since they are formed out of first (and possibly 
second) generation star ejecta. This may explain why the 
s-process enhanced stars do not show the same central 
concentration compared with the He-enriched stars, and 
why the s-process enhancements seem to be increasing 
already at relatively low [Fe/H], not much higher than 
that of the first generation . As a matter of fact, the 
data of Marino et al. (2012 Fig. 6 & 7) indeed shows two 
La-enrichment sequences, one belonging to what is iden- 
tified as the 'first generation' that correlates with C, and 
one that correlates rather with N, and associated with 
the second generation. In our scenario, the C-rich, low 
[Fc/H] and high La-abundan ce st ars would in fact be the 
fourth generation (see Sect. 5.2 for more discussion on 
this point). 

As discussed in Sect. T2J21 some w orkers (e.g. D'Ercole 
et al.|2008| |Conroy fc Spergel||2011[ ) favour the idea that 
the purging of gas from clusters is mostly due to super- 
nova. Massive-star SN II explosions would clear out the 
gas until a lower mass limit is reached where either SN 
lis no longer occur, or the explosions are so weak that 
they are incapable of removing all of the gas. Whether 
or not SN are effective in purging the gas sufficiently to 
enable the formation of a second generation out of the 
isolated ejecta of stars from a particular mass range will 
not only depend on the total mass of the system and 
the energy of the supernova, but also on their location 
within the cluster (progenitor). While SN purging clearly 
seems to be a viable mechanism in some (possibly less 
massive) clusters, our proposed scenario appears to be 
capable of explaining the detailed properties of the very 
unusual case of u> Cen, such as the small difference in 
[Fe/H] between the first and second stellar generations, 
and the larger difference in metallicity between the sec- 
ond and third generations. Moreover, it should also work 
for other clusters that are massive enough to retain SN 
ejecta. 

5. CONCLUSIONS 

5.1. Implications of the Galactic plane passage gas 
purging model 

As~ 
(like 



D'Ercole et al. (20081 put it, He-rich populations 
he bMS in u> Cen) are only the tip of the iceberg of 
the phenomenon of second (or third, etc.) stellar popula- 
tions. While our scenario of periodic gas purging events 
caused by the passages of clusters through the Galactic 
plane appears to be able to explain quite naturally the 
origin of the helium-rich bMS in uj Cen (and possibly 
the more Fe-rich MS-a/RGB-a population, as well), one 
can readily imagine that this simple picture, if applied to 
clusters having a wide range of progenitor histories (e.g., 
masses and orbits), may, in fact, generate quite a variety 
of realizations in the globular clusters that we observe 
today. Since each population would be associated with 
the ejecta arising from stars within one, or very few, rel- 
atively narrow initial mass ranges, isolated by successive 
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Table 5 

Summary of the f our generations ident ified in the Galactic plane passage gas purging model for u) Cen, and how they match with the 
cluster analysis of |Gratton et al~| ( |2011[ l and with the observational properties re porte d in DOrazi:2011jf and |Marino et aT] ( |2012^ (see 

text for details, in particular Sect.|5.3[l, 



generation parent expected abundance markers 
genera- 
tion 



Gratt onl co rresponds to feature s reported Leiden star numbe rs 
let al.l inlMarino et al.M2012|) ID'Orazi et al.N2011| 

j2011| 
group 

16015, 37247, 33011, 

low O, high Na (most extreme?) 41039, 46092, 34029 
for medium [Fe/H] 

highest Na & La for highest 60066, 60073, 34180, 

[Fe/H] 48323, 54022 

O, Na, N, C - La relation of O- 44462 
poor/Na-rich group 



first 


proto- 


low Fe; normal He, La; high a, high 


#4, #6 




cluster 


(from r-process) Eu 




second, 


first 


low 0, intermediate Fe; high He, N, Na, 


#3, 


bMS 




n-capture 


#2a 


third, MS- 


second 


high Fe, C+N+O; higher N, He, n- 


#2b 


a/RGB-a 




capture 




fourth 


first (sec- 


low (medium) Fe; normal He, medium 


#5, #1 




ond) 


O, Na, high C, n-capture 





Galactic plane passages that clear out any gas that had 
accumulated since the last passage, any two (or more) 
populations that are present in different clusters should 
not be identical, though possibly similar. The contribut- 
ing segments of the IMF are statistically selected accord- 
ing to the progenitor orbit and the merger history. 

However, in our scenario, Galactic plane passage inter- 
vals of 40 to 100 million years do favor the formation of a 
second, helium-rich pop ulation. If the delay is 40 million 
years (as also found by iMarino et al. 20121, upper and 
lower mass limits can both be imposed by the Galactic 
plane passages. For time intervals as large as 100 million 
years, only the lower mass cut-off for the progenitor stars 
of the second generation is set by the Galactic plane pas- 
sage, while the higher mass cut-off would be due to SN 
purging. In fact, this case is s imilar to the model sug- 



gested b y|D'Ercole et al.||2008| Sect. 3.1), who adopt 100 
Myr for their parameter Atf. They do not mention what 
sets this time interval, aside from noting that, if Atf was 
longer, SN la and C-producing AGB stars would con- 
tribute to the second-generation star formation (see the 
above discussion). This is indeed why Atf should not 
be longer if the desired outcome is a helium-rich second 
generation, but we suggest that the main reason is that 
the Galactic plane passages which clear out the gas and 
terminate star formation occur every ^ 100 Myr. 

In summary, the simple principle of periodic Galactic 
plane passage purging in combination with low-velocity 
winds from massive AGB stars AGB wind predictions 
and their preferential cooling properties, may be able to 
account for the bMS abundances including its central 
concentration, the MS-a/RGB-a abundance patterns as 
well as the homogeneous distribution of and wide range 
of [Fe/H] observed in s-enhanced stars. 

5.2. O-Na anti-correlation, n-capture element 
abundances, and other observed properties of ui Cen 

One of the (possibly defining, Carretta et al.|2010 1 fea- 
tures of globular clusters is the more or less complete 
presence of the O-Na anticorrelation. How does it fit in 
with the Galactic plane passage gas purging model? The 
key question is whether (a) the anti-correlation occurs 
within a sub-population, or (b) it is the superposition of 
the rather distinct O-Na abundances of the present sub- 
populations. The answer will depend on how the latter 
are i dentified. If one uses o nly the [Fe/H] abundance 
(e.g. D'Antona et al. 2011b) one may indeed combine 



N-rich/C-poor with C -rich/N-poor stars into one group 
(Marino et al. 2012 Fig. 2) and conclude that abun- 
dance anti-correlations are present even within a sub- 
population. 

This view point may motivate, or even require as an 
expla nationof the Na-0 anti -correlation, a dilution sce- 
nario ( |D'Ercole et"aL|2011b| ) which assumes that, while 
second-generation stars form in the cluster, unprocessed 
pristine gas is accreted from outside of the cluster. How- 
ever, if option (b) is in fact the case, for example, if the 
sub-population identification purely by [Fe/H] is not en- 
tirely accurate, then an alternative interpretation of the 
anti-correlation is possible in which dilution may not play 
an important role. 

Option (b) may also be favoured by observations of 
rather uniform Al abundances in the SGB-a popul ation 
in our scenario th e third generation, see Table[5]) by jPan" 



cino et al. (2011), suggesting that no anti-correlation is 
present m this sub-population. Although we have tried 
to resist the temptation in this paper to apply the Galac- 
tic plane pa ssage gas purging mo del to other clusters 
we note that |Carretta et aL (2012) found that the anti- 
correlations involving Al, Mg, IN a and O in NGC 6752 
- another GC for which multiple populations possibly 
including large He-enrichements have been found - man- 
ifest themselves in a rather discrete fashion, in which dif- 
ferent levels of enhancement and depletion cluster around 
discrete values that can be associated with individual 
sub-populations, rather than a continuous distribution 
that would be sugggested by a pure dilution mechanism. 

The Galactic plane passage gas purging scenario im- 
plies that two distinct populations may have the same 
[Fe/H] abundance, and therefore that a collection of stars 
with the same [Fe/H] may not necessarily belong to the 
same population. As a matter of fact, many possible 
elemental markers if used by themselves may lead to de- 
generate grouping of stars. In order to identify those 
stars that most likely represent a coevally formed sub- 
population, several elemental markers should be com- 
bined. Such a "group" analysis, taking into account 
four abund ance features simultan eously, has been per- 
formed by Gratton et al. (20111. The resultant sub- 
population identification is indeed rather suggestive of 
option (b) mentioned in the previous paragraph (super- 
position of the distinct O-Na abundance markers from 
sub-population). Of course, there is always a concern 
that the particular choice of group-finding criteria biases 
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the process . For that reas o n the sub-population identi- 
fication by Gratton et al. (2011) may be evolve in the 
future when more observational data are added to this 
kind of analysis. In any case, this approach seems to be 
an improvement over just using [Fe/H] to identify sub- 
populations. As a consequence, we have no compelling 
need for the dilution mechanism in our model for uj Cen. 

Our scenario has a few more implications - some of 
which are summarized in Tabled]- that we would like to 
briefly discuss. It offers an alternative interp retation of 
the star formation time-scale determination of lD'Antonal 



|et~al.| fgOllbJ that was based on the assumption that Fe 
in the MS-a/RGB-a populations originates in SN la. In 
our scenario, Fe for the third generation (MS-a/RGB-a) 
comes from the lowest mass, second genera tion S N II af- 
ter the second Galactic plane passage (Sect. 4~2~4| , which 
would be in agre ement with the a-ele ment abundance 
patterns found by 



rich population la 



Gratton et al.| (|2011|) for the most Fe- 
beled #2b. 



The time-scale limit from 
SN la would therefore not apply. In addition, star for- 
mation can take place even after the formation of this 
third generation from the late stellar winds of low-mass 
first-generation stars. 

In fact, if we retain our chemical evolution time 
range determination from CMD considerations (^ 1 Gyr, 
Sect. [2]), stars in the fourth generation can form from the 
slow wind ejecta of first-generation stars down to 1.8 M . 
For stars with < 2.6 M Q these ejecta would be C-rich 
(as well as La- rich from the main s process). Because 
the first-generation donor stars are not centrally concen- 
trated, the fourth generation stars that form out of them 
would also not be centrally located — if the cooling effi- 
ciency of lower-mass AGB ejecta are smaller, or because 
the mas s of the cluster is already sma ller at this point. 
Indeed, Johnson & Pilachowski (20101 report no radial 
gradient t or |La/Fe|. In any c ase, the observations re- 
ported by Marino et al. (2012) show that there are two 
different La-enrichment sequences, one associated with 
increasing C (O-rich/Na-poor group) and one associated 
with low O and high Na. In our scenario the former 
represents the transition from first to fourth generation, 
while the latter is associated with both the second and, 
possibly for even more extreme levels of enrichment, the 
third generation. 

If we interprete the c orrelation of C with La in the 
O-rich/Na-poor group of Marino et al. (2012) as the re- 
sult of the third dredge-up in 1.8 - 2.6 M first (and 
maybe second slightly higher [Fe/H]) generation stars, 
then their Fig. 7 implies that O and Na could also be 
the result of the third dredge-up. In fact, the slight in- 
crease of O with [Fe/H] in the O-rich/Na-poor group has 
been perceived to be quite a puzzle tha t implies an ex- 



tra source of O. However, the models of Herwig (2004b 
Fig. 6) for [Fe/H] = -2.3 do indeed predict [C/Fe] = 3.0, 
[O/Fe] = 1.8 and [Na/Fe] = 1.3 for the average abun- 
dance in ejecta of a 2 M Q model. The O comes from the 
dredge-up of primary He-burning products that becomes 
appreciable at these low metal abundances, and such O 
enhancements are indeed observed in many CEMP stars 
that may carry the mass transfer signature fr om gen- 



uine low-metallicity AGB stars (for example, Sivarani 



et al.| 2006| |Kennedy et al.||2011[ ). The Na in the 21V1 
model comes from s process in the He-shell flash convec- 
tion zone, where the neutron from the 22 Ne(a, n) 25 Mg 



reaction is captured again by 22 Ne. The predicted en- 
richment levels suggest that a population which forms 
out of such ejecta should have about twice the enrich- 
ment of C compar ed to O and Na, wh ich seems to be 
consistent with the Marino et al. (2012) data. 

Star Leiden 4446 2 has been tentatively identified by 
D'Orazi et al. (2011 ) as a mass transfer object in order to 
account for the extremely high C abundance. Although 
radial velocity measurements may support this possibil- 
ity, an alternative interpretation is that this star is part 
of the fourth generation (Table[5]) which forms out of the 
slow winds of 1.8 - 2.6 M Q first (or second for higher 
[Fc/H]) generation stars. In fact, this star coinci des very 
well with th e C-La correlation sequence shown in Marino 
et al.| ( 2012 Fig. 7, upper right pa nel). 



But the O-poor/Na-rich group of Marino et al. (|2012 
as wel l as all but the lowest [Fe/H] stars in D'Urazi et al. 
(2011 1, also show marked n-capture enhancements. Since 
in these stars the C abundance is low, the heavy elements 
cannot come from these lower-mass AGB stars. Instead 
we have to consider higher-mass AGB stars as well as 
super-AGB stars. T he n-capture element predictions of 



Karakas et al. (2012), e.g. their 6 M Q , Z = 0.0001 stellar 



evolution model, are based mostly on the 22 Ne neutron 
source in the Hc-shcll flash convection zone (although 
some contribution from a 13 C-pocket may be present 
as well). Clearly, the s-process models for these low- 
metallicity intermediate-m ass sta rs are quantiatively still 
ra ther uncertain ( c f. Sec t. 3.2.2 1. However, the models 
of Karakas et al. (2012) do predict that higher mass 
AGB stars at this metal content do produce s-process el- 
ements, possibly with significant enrichment factors, and 
with a ratio of light (Is: Sr, Y, Zr, Rb) to heavy (hs: Ba, 
La) s-process-elements that is higher than in the ejecta 
of lower mass stars. Such a signature would qualitatively 
agree with the n-capture abundances reported for the N- 
rich, intermediate- (bMS, second generation) a nd high- 
(MS-a/RG B-a, third generation) Fe group by (D'Orazi 
et al.||2011| . 

in addition to the n-capture production by 
22 Ne(a, n) 25 Mg in the He-intershell i n inte rmediate 
mass AGB stars, we mentioned in Sect. |3.3.l2| the pos- 
sible presence of i -process-conditions in super-AGB 
models with CBM (Herwig et al. 2012), which would 
provide for another non-standard source of n-capture 
elements in stars that produce He-rich and O-poor 
populations in uj Cen. Both of these sources would 
be responsible for the La-enhan ced sequence of th e 
O-poor/Na-rich group reported by Marino et al. (2012 1. 

An important consequence of this discussion is the no- 
tion that the Na enhancements which we expect from 
AGB stars from 1.8 M all the way up to the super-AGB 
stars at < 6.8 - 7.7 M Q may always be expected to go 
along with n-capture element enhancements. This assess- 
ment is supported by the La-Na correlations of both the 
O-rich /Na-p oor and the O-poor/Na-rich group of |Marino 
et al. (2012 lower-right panel of Fig. 7). 



5.3. Summary 

Obviously the observational identification of the four 
generations that we have specified - within the Galac- 
tic plane passage gas purging model - as the result of 
highly idealized processes, is complicated by interference 
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and superposition effects as well as contributions from, 
e.g., supernova purging, turbulent and tidal mixing, and 
cooling and mass loss flows, all of which are expected 
in a real cluster. In addition, observational uncertain- 
ties may cause some migration between observationally 
identified groups, which is in addition to the principal 
difficulty of deciding on the criteria and procedures that 
are use d to group stars (as discussed at the beginning 
of Sect. 5.2). These complications impose a limit to the 
accuracy that we can expect in how well observed prop- 
erties can be matched to the predictions of any scenario. 
The best that we can ask for at this point is rather qual- 
itative agreement, and with this goal in mind we have 
summarized the alignment of the Galactic plane passage 
gas purging model with some recently reported observa- 
tional properties of to Cen in Table[5] 

While most of the entries in Table[5| are based on our 
discussion in the previous sections there is a noteworthy 
peculiarity in the most metal-poor group ("first" gener- 
ation) which shoul d rep resent the genuine first genera- 
tion stars. In Sect. 15.21 we made the case that the ob- 
servational data may support the case of the superposi- 
tion of separate sub-populations forming the overall anti- 
correlation in uj Cen. It seems that, in fact, in this lowest 
[Fe/H] bin there are signs of an intrinisic anti-correlation 
with at most an unclear sign ature of n-capture en hance- 
ment. Groups #6 and #4 of Gratton et al. (2011 1 taken 
together display a significant scatt er in the [JNa/U ] ratio , 
while the lowest-metallicity bin of |D'Orazi et al.|(|2011|) 
includes two very N-rich stars (see also |Marino et al.| 
2012 1. This possible substructure in the most metal-poor 



population in uj Cen cannot be explained by the Galcatic 
plane passage gas purging model, and may be a relic of 
the proto-cluster object. 

This paper was significantly improved and extended as 
a result of the thorough and thoughtful report provided 
by the referee, Santi Cassisi, and we are very grateful to 
him for the time and effort that he put into his critique. 
We thank Antonio Sollima and Luigi Bedin for provid- 
ing the photometry of uj Cen that has been used in this 
investigation, as well as helpful comments. We are also 
grateful to Santi Cassisi for sending us his color trans- 
formation tables applicable to ACS photometry, and to 
Leo Michel-Dansac, who kindly assisted with the orbital 
integration of uj Cen. Thanks go, as well, to Kim Venn, 
David Hartwick, Pavel Denisenkov, and the entire "stars 
group" at the University of Victoria for many stimulating 
discussions concerning various topics of stellar evolution, 
globular clusters, and stellar populations in general. This 
work has been supported by the Natural Sciences and 
Engineering Research Council of Canada through Dis- 
covery Grants to FH, DAV and JFN. This research has 
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